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ABSTRACT 


The  purpose  of  the  experiments  and  calculations  reported  herein 
was  to  use  the  electron  beam  technique  for  rotational  and  vibrational  - 
temperature  and  density  measurements  in  the  free  stream  of  low 
density,  hypersonic  flows  of  nitrogen.  Beam  energy  was  varied  from 
10  to  20  kev  with  currents  of  1.  0  to  4.  0  ma.  Free-stream  rotational 
temperatures  were  in  the  range  from  30  to  200°K,  and  vibrational  tem¬ 
peratures  varied  from  300  to  3100°K,  and  free-stream  densities  were 
on  the  order  of  10~®  to  10“^  gm/cm^.  The  variation  in  the  flow  param¬ 
eters  was  obtained  through  the  use  of  four  different  nozzles;  each  was 
operated  at  two  different  conditions.  The  electron  beam  was  injected 
through  the  flow,  and  the  resulting  spontaneous  emission  of  light  was 
spectroscopically  analyzed.  Light  intensity  measurements  were  made 
photoelectrically,  and  signal  amplification  and  averaging  were  accom¬ 
plished  with  a  lock-in  amplifier.  Vibrational  band  intensities  were 
measured  by  electronic  integration  of  the  rotational  line  structure. 

The  results  of  the  rotational  temperature  measurements  are  in  good 
agreement  with  the  free-stream  temperature  values  predicted  on  the 
basis  of  conventional  aerodynamic  calibration  except  for  the  case  of 
the  Tunnel  L  10  N2  nozzle  for  which  there  is  an  approximately  24-percent 
discrepancey.  Measured  values  of  gas  density  are  within  15  percent  of 
the  theoretical  predictions,  and  the  results  of  vibrational  temperature 
determinations  are  consistently  on  the  order  of  10  percent  higher  than 
the  expected  values. 
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1.1  TUNNEL  DESCRIPTIONS 

Low  Density  Hypersonic  Wind  Tunnels  (L)  and  (M)  are  continuous 
flow,  low  density,  hypervelocity  wind  tunnels  (Figs.  1  and  2,  Appendix  I), 
and  thorough  descriptions  of  these  facilities  are  contained  in  Refs.  1 
and  2.  Briefly,  each  tunnel  consists  of  a  direct -current  electric  arc 
heater,  a  stilling  chamber,  an  expansion  nozzle,  a  test  section  with 
probe  traversing  mechanism,  a  vacuum  pumping  system,  and  necessary 
instrumentation.  Tunnel  L  may  also  be  operated  in  a  cold-flow  condi¬ 
tion,  i.  e. ,  without  the  arc  heater.  Normal  tunnel  operating  conditions 
of  the  two  tunnels  are  within  the  following  ranges: 


Parameter 

Engineering  Units  (U.  S. ) 

MKS  Units 

Working  fluid 

Nitrogen  (N2) 

— 

Total 

temperature 

540  to  7200°R 

300  to  4000°K 

Mach  number 

4  to  18 

— 

Stagnation 

enthalpy 

130  to  2140  Btu/lbm 

72  to  1189  kgm  cal/kgm 

Velocity 

2000  to  8500  ft/sec 

610  to  2590  m/s 

Static 

pressure 

0.  006  to  0.  095  torr 

— 

Static  density 

3  x  10"6  to  4  x  10"5  lbm/ft3 

4.  8  x  10"5  to  6.  4  x  10"4 

Unit  Reynolds 
number 

40  to  200  per  in. 

1575  to  7874/m 

Diameter  of 

uniform  core 

0.  5  to  10  in. 

0. 0125  to  0. 254  m 

Both  tunnels  are  facilities  for  investigating  the  simulation  of  flight 
at  extreme  altitudes  and  are  able  to  provide  aerodynamic  coefficients 
for  lifting  and  ballistic  bodies  in  the  hypersonic,  cold-wall  condition  at 
simulated  altitudes  above  300,  000  ft. 


1.2  GENERAL  EXPLANATION  OF  THE  ELECTRON  BEAM  TECHNIQUE 

A  method  for  directly  measuring  nitrogen  rotational  and  vibrational 
temperatures  and  number  density  in  wind  tunnel  flows  is  the  electron 
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beam  technique  pioneered  by  E.  P.  Muntz  (Refs.  3,  4,  5,  and  6).  This 
technique  involves  the  ionization  and  excitation  of  nitrogen  molecules 
with  a  beam  of  moderately  energetic  electrons  (10  to  100  kev).  The 
resulting  spontaneous  emission  in  the  spectral  region  from  3500  to 

o 

5000  A  at  low  nitrogen  density  is  predominately  that  of  the  first  nega¬ 
tive  system  of  nitrogen  (N2+  [l-]).  The  intensity  distribution  of  this 
system's  vibrational  bands  can  be  related  to  an  effective  vibrational 
temperature  of  nitrogen,  and  the  rotational  structure  of  the  vibrational 
bands  has  an  intensity  distribution  that  can  be  related  to  an  effective 
rotational  temperature  of  nitrogen.  Additionally,  the  intensity  of  any 
one  vibrational  band  can  be  related  to  the  density  of  nitrogen  molecules. 

At  low  gas  densities  the  beam  is  spatially  well  defined;  therefore, 
observation  of  emission  at  different  positions  along  the  beam  allows 
point  temperature  measurements. 

By  photoelectric,  spectroscopic  measurement  of  the  rotational  line 
intensities  of  a  vibrational  band,  a  rotational  temperature  can  be  deter¬ 
mined.  Vibrational  band  intensities  can  be  measured  by  electronic 
integration  of  the  rotational  line  intensities  of  a  vibrational  band. 

Usually  the  intensities  of  two  bands  are  used  to  form  a  ratio  which  can 
be  compared  to  a  theoretical  prediction  of  the  ratio  versus  vibrational 
temperature. 

For  the  density  range  of  interest  the  intensity  of  a  vibrational  band 
is  a  linear  function  of  density.  Therefore,  determination  of  the  intensity 
of  a  vibrational  band  at  some  unknown  density  relative  to  the  intensity  at 
a  condition  of  known  density  and  temperature  allows  the  gas  density  to  be 
determined  provided  the  vibrational  temperature  is  known. 


SECTION  II 

TUNNEL  INSTRUMENTATION  AND  OPERATING  PROCEDURE 


2.1  INSTRUMENTATION 

Figure  3  indicates  the  measurements  that  were  necessary  to  en¬ 
sure  that  proper  flow  conditions  were  obtained  in  the  tunnels.  Mass 
flow  rates  m^  and  m2  were  obtained  from  pressure  readings  from 
P^B  transducers.  The  input  temperature  of  the  nitrogen  flow  (TC^) 
was  obtained  from  a  copper -constantan  thermocouple  reading.  The 
stilling  chamber  pressure  (pQ)  was  read  using  a  P^B  transducer. 
Centerline  impact  pressure  (p')  was  monitored  with  a  ±0.  15  psid  trans¬ 
ducer  in  Tunnel  L  and  with  a  30-mm  Baratron®  pressure  gage  in 
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Tunnel  M.  Normally  the  stagnation  probe  was  at  a  position  of  1.5  in. 
from  the  nozzle  exit  in  Tunnel  L.  and  14.  25  in.  from  the  nozzle  exit  in 
Tunnel  M.  The  test  cabin  ambient  pressure  of  Tunnel  L  was  monitored 
with  the  ±0. 15  psid  transducer  and  with  an  Equibar®  pressure  gage;  that 
of  Tunnel  M  was  monitored  with  a  3-mm  Baratron  pressure  gage.  Arc 
heater  current,  I,  and  voltage,  E,  were  monitored  with  standard  meters. 
More  detailed  discussion  of  the  tunnel  instrumentation  may  be  found  in 
Ref.  2. 


2.2  OPERATING  PROCEDURE 

Before  beginning  each  tunnel  operation  the  pressure  transducer 
readouts  and  amplifiers  were  balanced,  and  for  Tunnel  L  the  ±0.  15  psid 
transducer  was  calibrated  over  the  range  of  expected  p^  and  prp  values 
using  a  micromanometer  to  read  the  calibration  pressure.  Mass  flow 
rates  m^  and  m2  were  then  set  to  the  required  values  and  pQ,  TC^,  m^, 
and  m2  were  read. 

For  arc -heated  flow  the  heater  was  then  started,  and  the  power 
setting  adjusted  to  obtain  the  required  pQ.  Values  of  E,  I,  p^,  and  pT 
were  then  read  to  ensure  proper  flow  conditions. 


SECTION  III 

ELECTRON  BEAM  SYSTEMS 


3.1  VACUUM  SYSTEM 

In  order  to  inject  the  electron  beam,  emitted  by  the  electron  gun  at 
a  pressure  of  less  than  1.  0  x  10"^  torr,  into  the  higher  free-stream 
pressures  of  both  tunnels,  a  dynamic  pumping  system  was  used. 
Schematic  diagrams  of  the  system  are  shown  in  Figs.  4,  5,  and  6.  The 
system  consisted  of  a  4 -in.  oil  diffusion  pump  with  a  hlank-off  speed  of 
approximately  700  liters /sec  and  a  mechanical  forepump  with  a  maxi¬ 
mum  pumping  speed  of  140  liters /min.  An  additional  mechanical  pump 
with  70-liters /min  pumping  speed  was  used  for  rough-pumping  the  sys¬ 
tem.  The  electron  gun  section  was  separated  by  an  orifice  from  the 
tunnel  free  stream  into  which  the  beam  was  injected.  The  orifice  diam¬ 
eter  was  nominally  0.  040  in.  and  the  orifice  length  was  0.  25  in.  Pres¬ 
sures  in  the  electron  gun  section  were  normally  maintained  by  this  sys¬ 
tem  at  less  than  2  x  10-5  torr  as  measured  with  an  ion  gage. 
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As  shown  in  Figs.  5  and  6  the  electron  gun  was  mounted  in  a  housing 
inside  the  tunnel  test  cabin.  Atmospheric  pressure  was  always  main¬ 
tained  inside  the  gun  housing  to  prevent  high  voltage  arcing.  The  gun 
housing  was  bolted  to  a  traversing  mechanism  that  allowed  axial  position¬ 
ing  of  the  electron  beam  along  the  tunnel  free  stream.  In  order  to  pre¬ 
vent  accidental  overpressuring  of  the  electron  gun  which  would  shorten 
considerably  the  gun  lifetime,  an  electrically  operated,  pneumatic  gate 
valve  was  included  between  the  orifice  and  the  electron  gun. 


3.2  ELECTRON  GUN  AND  RECEIVER  CUP 

The  electron  guns  were  television-type,  oxide  cathode  models  modi¬ 
fied  for  a  maximum  of  50-kv  operation.  The  accelerator  grid  of  these 
guns  was  also  modified  by  having  its  opening  enlarged  to  3/ 16-in.  diam¬ 
eter.  Magnetic  focusing  and  deflection  were  used  and  the  filament  was 
battery  powered  at  6  to  18  v. 

The  electron  gun  high  voltage  power  supply  provided  0  to  50  kv, 

0.  01 -percent  ripple  voltage,  at  0  to  5  ma.  The  high  voltage  supplied  to 
the  guns  was  always  negative  with  respect  to  the  anode  which  was  at 
ground  potential. 

The  electron  gun  system  was  capable  of  injecting  up  to  a  4.  5-ma 
beam  at  10  to  30  kv  across  the  tunnel  free  stream  to  a  beam  receiver. 
The  electron  beam  receiver  cup  for  Tunnel  L  is  shown  in  Fig.  7. 
Basically  it  was  a  12-in.  length  of  1-3/4-in,  -diam,  water-cooled, 
stainless  steel  tubing  capped  with  a  stainless  steel  plate.  The  double 
grid  structure  added  another  inch  to  the  cup.  The  grids  were  1/16 -in. 
mesh  made  of  0.  02-in.  -diam  brass  wire.  The  entrance  diameter  of 
both  grids  was  1  in. 

J 

The  interior,  or  second,  grid  was  negatively  biased  at  25  to  30  v 
to  prevent  the  ejected  secondary  electrons  from  leaving  the  cup  interior, 
which  was  carbon  coated.  The  outer  grid  and  cup  are  electrically  con¬ 
nected,  and  the  entire  cup  was  electrically  isolated  from  the  tunnel 
floor  by  two  sheets  of  Mylar®  film. 

Beam  current  measurement  in  Tunnel  M  presented  more  of  a 
problem.  A  receiver  cup  of  the  type  used  in  Tunnel  L  would  have  its 
entrance  3-1/2  ft  from  the  orifice,  which  presented  the  problem  of 
collecting  the  entire  beam.  For  vibrational  and  rotational  temperature 
measurements  a  precise  value  of  the  current  is  unimportant,  and,  for 
all  temperature  data,  was  a  nominal  value  as  registered  by  an  aluminum 
flat-plate  collector  of  approximate  dimensions  12  by  18  in.  located  on 
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the  tunnel  floor,  approximately  4  ft  from  the  orifice.  The  accuracy  of 
such  a  measurement  is  indeed  questionable  because  of  both  boundary 
layer  and  secondary  ejected  electron  effects;  however,  the  value  of 
current  was  constantly  monitored  to  ensure  minimum  current  fluctua¬ 
tions  and  drift,  thereby  decreasing  the  imprecision  of  the  temperature 
measurements.  As  shown  later,  knowledge  of  the  value  of  the  beam 
current  is  critical  for  an  accurate  density  measurement,  and  for  that 
reason  it  was  deemed  essential  that  an  additional  current  measurement 
be  made  in  Tunnel  M,  specifically  at  the  point  of  injection  of  the  beam 
into  the  flow,  thereby  minimizing  the  errors  attributable  to  beam 
spreading.  To  do  this  a  collector  was  mounted  onto  a  traversing 
mechanism.  Before  obtaining  an  absolute  density  data  point,  this  col¬ 
lector  was  traversed  axially  until  it  was  located  directly  beneath  the 
gun  orifice  (Fig.  6).  The  current  at  the  injection  point  was  measured, 
and  the  auxiliary  collector  was  withdrawn  to  allow  the  electron  beam 
to  traverse  the  tunnel  flow.  Figure  8  is  a  photograph  of  the  auxiliary 
collector  which  was  constructed  of  a  2 -in.  copper  elbow  with  a  re¬ 
flecting  plane  attached  to  give  an  effective  increase  in  the  length-to- 
diameter  ratio,  thereby  minimizing  the  loss  of  ejected  secondary 
electrons.  The  location  of  this  current  measurement  point  was 
approximately  24  in.  above  the  tunnel  centerline.  The  individual  values 
of  current  from  each  of  the  two  collectors  in  Tunnel  M  were  read  from 
a  digital  voltmeter. 


3.3  SPECTROSCOPIC  INSTRUMENTATION 

Schematic  drawings  of  the  spectrometer  setups  for  Tunnels  L 
and  M  are  shown  in  Figs.  9  and  10,  respectively.  In  the  Tunnel  L  set¬ 
up  a  14-in.  -focal  length  lens  focused  the  optical  radiation  produced  as 
a  result  of  electron  beam  excitation  onto  the  spectrometer  slit.  The 
f/numbers  of  the  lens  and  spectrometer  were  matched  at  f/8.  7.  When 
band  integration  measurements  were  made,  a  dove  prism  was  placed 
between  the  lens  and  spectrometer  slit  to  rotate  the  beam  image  per¬ 
pendicular  to  the  entrance  slit.  A  chopping  disk  was  located  near  the 
entrance  slit  for  modulating  the  light  input  at  400  Hz.  For  the  Tunnel  M 
setup  the  optical  radiation  was  gathered  by  an  f/12,  48-in.  -focal  length 
lens,  and  the  chopping  disk  operated  at  340  Hz. 

The  spectrometer  for  both  setups  was  a  1-meter  Jarrell-Ash 
Czerny-Turner  combination  scanning  spectrometer-spectrograph.  The 
102-  by  102-mm  diffraction  grating  used  had  1180  grooves /mm  and  was 
blazed  at  5000  A.  With  a  20-micron  slit  width  the  resolution  in  the  first 

O 

order  was  0.  29  A,  which,  although  by  no  means  the  maximum  attainable 
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resolution,  was  sufficient  for  resolution  of  the  R-branch  lines  for  most 
of  the  prominent  N24  [l“]  bands.  Scanning  speeds  were  2.  5  A /min  for 
rotational  spectra  and  25  A /min  for  vibrational  band  integration. 

A  photomultiplier  tube  mounted  at  the  exit  slit  of  the  spectrometer 
was  used  as  the  detection  element.  The  photo-cathode  voltage  applied 
to  the  tube  was  normally  -1700  vdc,  and  a  load  resistance  of  500  kilohms 
was  used  on  the  photomultiplier  output;  output  signals  were  processed  by 
a  lock-in  amplifier.  For  the  Tunnel  L.  system  the  reference  signal  to 
the  lock-in  amplifier  was  provided  by  a  small  6-v  light  bulb  and  photo¬ 
multiplier  tube.  For  the  Tunnel  M  system  the  reference  signal  to  the 
lock-in  amplifier  was  supplied  by  the  340-Hz  chopping  wheel  assembly. 

In  both  tunnels  the  output  of  the  lock-in  amplifier  was  recorded  by 
a  strip-chart  recorder  for  rotational  spectra  scans.  In  Tunnel  L  for 
band  intensity  measurements  the  output  of  the  lock-in  amplifier  was 
fed  to  an  active  electronic  integrator  (Fig.  11)  with  an  accuracy  of 
±1  percent  provided  the  integration  time  was  no  longer  than  about 
200  sec.  The  integrator  output  was  recorded  by  the  strip-chart  re¬ 
corder.  In  Tunnel  M  the  lock-in  amplifier  output  was  fed  to  a  different 
integrator  (Fig.  12)  with  the  output  displayed  on  a  digital  voltmeter. 
Because  the  amplifier  of  the  Tunnel  M  integrator  had  a  factor  of  300 
increase  in  gain  and  a  factor  of  100  drift  decrease,  the  Tunnel  M  inte¬ 
grator  was  more  accurate  than  the  one  used  in  Tunnel  L. 

The  relative  sensitivity  of  the  optical  systems  as  a  function  of  wave¬ 
length  was  measured  using  a  calibrated,  tungsten-strip  standard  lamp. 
Sensitivity  curves  for  the  optical  systems  are  shown  in  Figs.  13  and  14. 

Because  of  the  small  uniform  flow  core  diameter  in  Tunnel  L, 
radial  scanning  of  the  flow  field  was  not  performed,  and  all  measure¬ 
ments  were  made  on  the  free -stream  centerline.  Optical  alignment  was 
accomplished  by  positioning  a  small  light  source  on  the  tunnel  centerline 
at  the  desired  axial  position  and  focusing  the  light  source  image  onto  the 
spectrometer  slit.  In  Tunnel  M  radial  scanning  was  performed  by 
mounting  the  lens  and  spectrometer  on  a  welding  table  which  was  capable 
of  being  raised  or  lowered,  thus  moving  the  spectrometer  and  optics 
system  along  the  radial  direction  of  the  beam.  To  ensure  that  the  elec¬ 
tron  beam  intersected  the  tunnel  centerline,  an  object  was  located  on 
the  centerline  and  focused  into  view  using  a  leveling  telescope  on  the 
side  of  the  tunnel  opposite  to  the  spectrometer  (Fig.  10).  A  plumb  bob 
was  then  suspended  from  the  electron  beam  orifice,  and  the  beam  mount 
was  rotated  until  the  plumb  bob  intersected  the  centerline,  which  was 
indicated  by  the  image  in  the  level.  The  lens  was  then  adjusted  until 
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the  entrance  slit  of  the  spectrometer  came  into  focus  with  the  level,  and 
at  that  point  it  was  certain  that  ( 1)  neglecting  magnetic  field  effects 
which  were  small,  the  electron  beam  intersected  the  centerline,  and 
(2)  the  spectrometer  entrance  slit  was  focused  onto  the  beam.  The 
radial  position  of  observation  was  determined  by  mounting  a  scale  on 
the  spectrometer  table;  using  the  level,  a  centerline  reference  position 
was  determined  on  the  scale,  which  had  a  1-mm  least  count,  and  the 
table  was  raised  until,  as  observed  through  the  level,  a  change  of  posi¬ 
tion  of  predetermined  magnitude  had  occurred.  After  the  magnetic  field 
coil  was  turned  on  during  tunnel  operation,  a  final  optimizing  adjust¬ 
ment  of  the  axial  and  angular  position  of  the  beam  was  required,  but  the 
change  was  always  small  and  did  not  require  refocusing;  both  tunnels 
required  this  adjustment. 


SECTION  IV 

ELECTRON  BEAM  TECHNIQUE 


4.1  INTENSITY  EQUATIONS 


On  the  basis  of  Muntz's  excitation  theory  the  intensity  equation  for 
the  rotational  lines  in  the  R  branch  of  a  N2+  [l“]  system  vibrational 
band  excited  by  an  electron  beam  may  be  written  as 


fn 


flv'.K'  ^ K  ' 

v  2,  K  -1  v  2.  K  -1 
K'(G>„ 


=  -B0  K'(K'+  1) 


he 


+  constant 


(1) 


in  which 


(G)„ 


(K'+  l)e-2B0K'(K'+l)hcATR  +  (K')e+2B0K 'hc/kTR 
2K'  1 


(2) 


Collisional  deexcitation  effects,  on  rotational  line  intensity  have  been 
neglected  in  the  development  of  Eq.  ( 1).  Furthermore,  it  has  been 
assumed  that  the  vibrational  temperature,  Tv,  is  sufficiently  low  to 
allow  the  neglect  of  all  ground  electronic  state  vibrational  levels  except 
the  lowest  (v^*  =  0). 

For  the  density  range  of  both  tunnels  the  first  assumption  should  be 
valid,  but  the  second  assumption  is  not  strictly  valid  under  the  high 
vibrational  temperature  conditions  of  both  tunnels.  However,  calcula¬ 
tions  have  shown  that  even  for  Tv  =  3100°K  the  measured  rotational 
temperature  differs  by  only  -2  to  -3°K  over  the  Tj^  range  from 
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50  to  300°K  from  those  values  of  Tj^  obtained  assuming  Tv  =  500°K  or 
less,  thereby  justifying  the  neglect  of  Tv  effects  on  Tj^  measurements. 


As  shown  in  Ref.  7  the  intensity  of  a  vibrational  band  is 


I(v  V  2'  )  = 


£  p<V.M'>e-G(T1'',hc/kT« 

(const.)  Nt  p(v  v  2')  »  j  "■=  0 


(p.M, 


X4  ' 


.  v2 


OO  /  /  //v 

m  P(v  .  V2  ) 

r2*0 


Then  the  intensity  ratio  of  two  bands  (a^,  bv2*>  and  (cv> ,  dv2*)  is 


(3) 


I(a,',br'') 

Hcv',dv^) 


-  /  ^cd V  P<av'  by;*)  t{ 
\^«b/  p(cv*  dvT)  1 


=  0 


£  p<cv',  yl")e“G<Tl">hc/l‘TV 

//  _ 

VI  =0 


—  OO 

— 

£ 

**  _ 

p(c¥',  v2")/A.3Cv^  „2" 

v  2  =0 

OO 

£ 

P<«v  V2')/A3Bv',  v2' 

v2  =0 

__ 

(4) 


Here  again  collisional  deexcitation  effects  have  been  neglected. 


4.2  ROTATIONAL  TEMPERATURE  MEASUREMENT 


From  a  rotational  spectral  scan  the  peak  heights  of  the  spectral 
lines  are  measured  and  used  as  the  relative  intensity  values  of  the 
v'  K' 

rotational  lines,  I  ’  x,,  :  using  these  values  and  the  known  values 

v2*,  K  -1' 

of  K'  »  may  be  solved  for  rotational  temperature  by  an 

v2’*,  K'-l 


iterative  machine  calculation.  The  program  (Program  TROT)  for  the 
computer  is  shown  in  Table  I,  Appendix  II.  The  data  deck  requires  in¬ 
put  data  cards  for  the  following:  spectral  scan  number  and  designation, 
total  number  of  spectral  lines  in  a  scan,  and  the  relative  intensity  value 
of  each  line.  The  intensities  of  even-numbered  rotational  lines 
(K'  =  2,  4,  . . . )  are  doubled  before  input  to  account  for  the  influence  of 
nuclear  spin  statistics. 
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Letting 


x  =  K'(K'  +  1)  (6) 

the  machine  calculation  performs  a  least -squares  fit  of  the  y  values  to 
the  straight  line  given  by  Eq.  (1).  As  shown  from  Eq.  (1)  the  slope  of 
this  straight  line  is  hcB0/kTj^;  hence,  a  rotational  temperature  can  be 
calculated  from  the  slope  of  the  straight  line,  A  rotational  tempera¬ 
ture  of  200°K  is  assumed  for  the  calculation  of  (G)Q  in  the  first  step  of 
the  iteration  process,  and  a  straight  line  is  then  fitted  to  the  y  values,  the 
slope  of  which  yields  a  new  value  of  the  rotational  temperature.  If  this 
calculated  value  of  rotational  temperature  is  not  within  0.  01  percent  of 
the  assumed  value  of  rotational  temperature,  the  calculated  value  is 
used  as  the  assumed  value  of  rotational  temperature  for  the  calculation 
of  a  new  set  of  (G)0  values,  and  the  least-squares  fit  of  the 'intensity  data 
according  to  Eq.  { 1)  is  repeated  until  the  prescribed  iteration  accuracy 
of  0.  01  percent  is  achieved. 

As  observed  by  Williams  (Ref.  7)  and  Ashkenas  (Ref.  8)  only  for 
temperatures  near  200°K  (pure  nitrogen  gas,  equilibrium  conditions)  is 
the  plot  of  y  values  versus  x  =  K'(K/  +  1)  a  good  straight  line;  that  is, 
only  with  these  restrictions  is  the  temperature  determined  from  the 
slope  of  the  line  not  a  function  of  the  number  of  rotational  lines  used. 

For  this  reason  a  temperature  Tpj(Kmax)  defined  by  the  maximum  num¬ 
ber  of  spectral  lines  (Kmax)  is  first  calculated  by  the  iterative  pro¬ 
cedure  described  in  the  preceding  paragraph.  Then  the  entire  process 
is  repeated  for  -1  tines,-  and  this  gives  another  rotational  tem¬ 

perature  value  Tp^K^^-l).  This  process  is  repeated  i  times  through 
Kmax  “i  +  l  =  Only  if  the  variation  of  Tr  with  the  number  of  rota¬ 
tional  lines  used  is  relatively  small  can  a  single,  well-defined  tempera¬ 
ture  be  obtained.  The  standard  deviation  of  the  temperature  is  printed 
in  °K  for  each  (Kmax  -i+1)  data  set.  This  value  is  computed  using  for¬ 
mulae  from  Beers  (Ref.  9). 

For  a  given  tunnel  condition  four  to  six  spectral  scans  are  taken  to 
comprise  one  temperature  determination,  and  these  scans  are  averaged 
within  Program  TROT.  The  output  of  this  averaging  process  is  pre¬ 
sented  as  an  average  ratio  of  rotational  temperature  to  predicted  free- 
stream  temperature  (Tr/T0)  and  average  standard  deviation,  ct,  for  a 
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given  value  of  the  number  of  spectral  lines  used  in  calculating  a  Tr  in 
a  fixed  set  of  data.  Each  individual  value  of  used  in  the  calcu¬ 

lation  of  the  average  value  of  Tr/T,,,  is  weighted  according  to  its 
standard  deviation,  o,  as  given  by 


which  is  by  definition  a  weighted  average  (Ref.  9),  where 

<J{  =  the  standard  deviation  of  the  ith  spectral  scan 

yi  =  determined  rotational  temperature  divided 

1  by  the  free-stream  temperature  for  the  ith 

spectral  scan  (8) 

s  =  the  number  of  spectral  scans  for  a  given  tunnel 
condition  (normally  4  to  6). 


By  definition  (Ref.  9)  the  standard  deviation  of  the  average  value  of  a 
set  of  experimental  data  is 


where 


(9) 


(10) 


It  has  been  observed  (Refs.  7  and  8)  that  the  degree  of  nonlinearity  of 
a  Boltzmann  plot  of  rotational  line  intensities  and  also  the  discrepancy 
between  measured  rotational  temperature  and  known  temperatures  in 
laboratory  conditions  is  a  function  of  the  density  and  temperature  of  the 
gas  into  which  the  electron  beam  is  injected.  Therefore,  it  is  neces¬ 
sary  to  empirically  correct  the  rotational  temperature  values  deter¬ 
mined  from  Program  TROT,  the  magnitude  of  the  empirical  correction 
being  density  and  temperature  dependent,  and  this  is  accomplished 
through  knowledge  of  the  approximate  density  and  temperature  of  the 
tunnel  flow  and  from  data  from  laboratory  experiments. 


4.3  VIBRATIONAL  TEMPERATURE  AND  DENSITY  MEASUREMENT 

As  shown  in  Ref.  10  the  proper  equation  for  the  band  intensity 
ratio  of  two  vibrational  bands  neglecting  quenching  effects  is  Eq.  (.4) . 
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As  pointed  out  in  Ref.  10  the  last  factor  of  Eq.  (4)  is  identically  equal 
to  the  ratio  of  the  lifetimes  of  the  av»  and  cv/  vibrational  levels.  Hence, 
Eq.  (4)  may  be  written 


Kav  bv2  1  /^cd  \ 

I(cv',  d„2")  \Aab/ 


£  p(av'.  v1'VC(vl",l’c/lTv 

vf=0 


£  p(cv'  vl")e~G(vl")hc/ltT» 
=0 


(ll) 


in  which 

and 


p(v',  V2  ")  =  q  (v'f  V2  ")  R|[fv',  v2") 


p(v',  vi  ")  =  q  (v'vi")  R|  [fv',  Vl"l 


(12) 

(13) 


Calculations  by  Hornkohl  have  yielded  values  for  q(v',  )  and 

q(v',  vj*  ),  and  experimental  data  and  calculations  by  Williams,  Lewis, 
and  Hornkohl  have  determined  values  of  p(v',  V2'').  These  parameters 
are  listed  in  Tables  II,  III,  and  IV.  The  same  experimental  data  also 
yielded  values  for  Ty,  _g,  Tv,  =  i>  TV'  =  2  6.58  x  10”®,  6.49  x  10”®,  and 

6.43  x  10”8  sec,  respectively.  Since  the  values  of  p(v',  vg"  )  are 
experimentally  determined  with  concomitant  errors  it  is  considered 
prudent  to  minimize  such  errors  by  taking  advantage  of  the  fact  that 
fav,,  bV2*  =  ?cv,,  dV2*  if  a  -  b  =  c  -  d,  thereby  yielding 

Re^rcd^  =  Re^rab)»  so  that  the  ratio  p(av/,  bV2*  )/p(cy, ,  dV2„  )  is  deter¬ 
mined  predominantly  by  the  ratio  of  calculable  parameters 
q(av/,  bV2*.)/q(cv* ,  dV2")-  Unfortunately,  no  work  has  been  performed 

to  determine  excitation  band  strength  parameters  p(v',  vj*  ),  and  one 
must  resort  to  the  usage  of  the  set  of  excitation  Franck-Condon  factors 
q(v',  vr» ),  thus  producing  a  possible  inaccuracy  in  using  Eq.  (11). 


The  intensity  of  vibrational  band  pairs  is  determined  by  elec¬ 
tronically  integrating  the  rotational  structure  of  the  vibrational  bands; 
by  comparing  the  experimental  band  intensity  ratio  to  a  plot  of  calcu¬ 
lated  band  intensity  ratio  versus  vibrational  temperature,  a  vibrational 
temperature  is  determined.  Plots  of  the  most  useful  band  intensity 
ratios  versus  Tv  are  given  in  Figs.  15  and  16. 

It  has  been  observed  in  Ref.  11  that  the  relationship  of  N2+  [l”l 
vibrational  band  relative  intensity  to  nitrogen  density  (Eq.  (3))  is  linear 
up  to  a  density  of  approximately  7.  5  x  10" ?  gm/cm^  which  includes  the 
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density  range  of  both  Tunnels  L.  and  M.  However,  it  can  be  seen  from 
Eq.  (3)  that  the  relative  intensity  is  also  a  strong  function  of  vibrational 
temperature.  Figures  17  and  18  show  the  relative  intensity  of  a  (0,  i) 
band  as  a  function  of  nitrogen  density  in  terms  of  equivalent  room  tem¬ 
perature  pressure  (ERTP)  and  vibrational  temperature.  The  nitrogen 
density  measurements  can  therefore  be  made  by  electronically  inte¬ 
grating  a  vibrational  band  intensity  at  the  desired  tunnel  condition  and 
comparing  this  to  the  integrated  value  at  a  known  condition  of  pressure 
and  temperature  and  correcting  for  vibrational  temperature  effects 
using  the  measured  value  of  vibrational  temperature. 


SECTION  V 

EXPERIMENTS  AND  RESULTS 


5.1  4  N2  NOZZLE-TUNNEL  L 

The  4  N2  nozzle  is  a  cold  flow  nozzle  operated  at  either  of  two 
fixed  conditions,  the  flow  parameters  of  which  are  listed  below: 


Mn  « 

Standard  Condition  -  -  -- 

p  =0.  0160  atm 
o 

T.  =  67°K  (120.  6°R) 

U, ,,  =  678.  2  m/sec  (2225  ft/sec) 
Pm  =  4.  8  x  IQ'7  gm/cm3 

(  3.  0  x  IQ'5  lbm/ft3) 
p  =  0.  075  torr 
=  4.  05 

ERTP  =  0.  320  torr 


.  Off -Standard  Condition 
pQ  =  0.  0184  atm 
T,,  =  66°K  (118.  8°R) 

Ub,  =  679.  7  m/sec  (2230  ft/sec) 
pa  -  5.  6  x  ID'7  gm/cm^ 

(3.5  x  10"5  lbm/ft3) 
p  =  0.  092  torr 
M*  =  4.  10 
ERTP  =  0.  371  torr 


The  foregoing  conditions,  and  those  presented  as  predicted  or  cali¬ 
brated  for  other  flows,  are  based  on  nozzle  calibrations  made  with 
impact  pressure  probes,  directly  measured  reservoir  pressure,  and 
total  temperature  derived  both  from  nonequilibrium  thermo-chemical- 
kinetic  calculations  and  from  measured  total  enthalpy  of  the  test  section 
flow.  With  and  pQ  determined,  the  ratio  p'  /pQ  yields  a  Mach  num¬ 
ber,  M,,,  for  a  particular  gas.  Then  T,,,,  p^,  and  pa  are  calculated  on 
the  basis  of  pQ,  TQ,  and  M,,,  using  the  same  gas  thermo-chemical-  kinetic 
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model  assumed  in  deriving  M0  from  /pQ.  This  model  is  a  perfect 
gas  Tor  the  MB  =  4  case.  Thus,  T,,,  and  pa  are  not  directly  measured  in 
the  gas -dynamic  flow  calibration,  but  their  values  should  be  reasonably 
certain  at  M,,,  =  4  because  of  the  accuracy' of  p^ ,  pQ,  and  T0  values  and 

the  rather  strong  basis  for  the  assumption  of  a  frozen  (7  =  1.  400)  expan¬ 
sion.  For  the  nozzle  flows  to  be  discussed  later,  the  same  could  be 
said,  except  then  it  is  possible  that  the  flow  does  not  freeze  to  a  constant 
7  process  until  it  reaches  the  nozzle  throat.  The  freezing  referred  to 
concerns  molecular  vibration  which  is  believed  the  only  nonequilibrium 
mode  in  any  case  considered  (Refs.  1  and  12). 

If  appreciable  error  exists,  in  the  higher  TQ  flows,  the  worst  cases 
studied  show  that  the  values  of  T0  predicted  from  the  standard  nozzle 
calibrations  may  be  low  by  as  much  as  30  or  40  percent,  whereas  the 
predicted  values  of  p^  should  be  high  by  no  more  than  5  or  10  percent. 
These  approximate  maximum  uncertainty  levels  are  based  on  com¬ 
parison  of  the  listed  Tm  and  p a  values  with  corresponding  results  which 
would  have  been  calculated  if  full  thermo-chemical-kinetic  equilibrium 
had  been  assumed  throughout  the  nozzles  (cf. ,  Ref.  l). 

Beam  measurements  were  made  approximately  1  in.  from  the  noz¬ 
zle  exit  at  =  4.  The  length  of  the  beam  observed  was  approximately 
0.  25  in.  in  the  radial  direction  centered  on  the  flow  core  centerline, 
and  the  diameter  of  the  flow  core  was  approximately  1.  25  in. 

Two  sets  of  six  rotational  line  scans  were  made  for  each  nozzle 
condition.  A  typical  spectral  scan  is  shown  in  Fig.  19,  and  the  tem¬ 
perature  results  (calculated  by  Program  TROT)  are  shown  in  Figs.  20 
and  21.  The  variation  of  the  measured  Tr  value  with  rotational  quantum 
number  and  the  discrepancy  between  the  measured  value  and  the  pre¬ 
dicted  value  were  to  be  expected  on  the  basis  of  the  results  in  Refs.  7 
and  8  (see  Figs.  22  and  23).  From  Fig.  23  it  is  seen  that  a  13  to  20°K 
correction  to  the  experimental  data  is  needed  when  ten  rotational  lines 
are  used  in  the  Tjj  determination.  From  Figs.  20  and  21  this  would 
yield  Tjj  values  from  60  to  67°K.  For  the  case  of  15  rotational  lines 
used  to  determine  Tr  a  28  to  32°K  correction  must  be  made.  From 
Figs.  20  and  21  this  would  yield  values  from  64  to  68°K.  Therefore, 
reasonable  agreement  between  measured  and  predicted  Tp  values  was 
obtained  from  the  data  by  taking  into  account  effects  of  density  and  num¬ 
ber  of  rotational  lines  used. 

The  nitrogen  density  measurements  were  made  by  electronically 
integrating  the  N24  [l“]  (0,  0)  band  intensity  at  the  tunnel  condition  and 
comparing  this  to  the  integrated  value  at  a  static  condition  of  0.  100  torr 
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and  room  temperature.  The  measured  pm  values  were  5.  61  x  10"7  gm/cm3 
and  6.  56  x  10“ 7  gm/cm3  for  the  standard  and  off-standard  conditions, 
respectively.  The  measured  values  are  observed  to  be  approximately 
15  percent  high,  and  this  discrepancy  may  be  the  result  of  beam  spread¬ 
ing  effects  at  the  collector  cup  since  a  beam  of  only  20  kv  was  used.  It 
is  felt  that  a  repetition  of  those  measurements  using  a  beam  energy 
from  50  to  100  kev  would  prevent  significant  beam  spreading  and  might 
result  in  a  smaller  discrepancy  in  the  density  value. 

For  each  nozzle  flow  a  scan  of  the  beam  induced  radiation  was  made 
from  3500  to  6000  A  to  detect  the  presence  of  contaminative  radiation 
which  might  affect  the  measurements.  The  4  N2  nozzle  beam  radiation 
scan  is  shown  in  Fig.  24,  and  it  is  shown  that  all  significant  radiation  is 
attributable  to  excited  or  ionized  nitrogen. 


5.2  9  N2  NOZZLE-TUNNEL  L 

The  9  N2  nozzle,  designed  for  arc-heated  flow,  is  operated  at  either 
of  two  fixed  conditions.  The  flow  parameters  which  are  listed  below  are 
also  based  on  the  calibration  with  aerodynamic  probes  as  described  in 
the  previous  section. 

Mro  -  9 


Standard  Condition 

Off -Standard  Condition 

pQ  =  2.  04  atm 

Po 

=  1.  70  atm 

Tq  =  2365°K  (4257°R) 

To 

=  1662°K  ( 299 2°R) 

T*.  =  135°K  ( 243°R) 

Tm 

=  93°K  ( 167°R) 

U,,,  =  2204*m/sec  (7230  ft/sec) 

Uc 

=  1844  m/sec  (6050  ft/sec) 

pro  =  2.  02  x  10~7  gm/cm3 

Pco 

=  2.  32  x  10  ^  gm/cm3 

(126  x  10"5  lbm/ft3) 

(1.  45  x  10"5  lbm/ft3) 

p  =0.  0615  torr 

p. 

=  0.  049  torr 

M.  =  9.  37 

M_ 

GO 

=  9.  30 

ERTP  =  0.  134  torr 

ERTP 

-  0.  155  torr 

Again,  beam  measurements  were  made  approximately  1  in.  from  the 
nozzle  exit.  The  diameter  of  the  flow  core  was  approximately  2.  0  in. 

and  a  0.  25 -in.  beam  length  at  the  core  centerline  was  observed. 

Two  sets  of  six  rotational  scans  were  made  for  each  nozzle  condi¬ 
tion.  Typical  spectral  scans  are  shown  in  Figs.  25  and  26.  The  results 
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(calculated  by  Program  TROT)  are  shown  in  Figs.  27  and  28.  On  the 
basis  of  the  laboratory  data  findings  shown  in  Fig.  23  for  78°K  data  and 
Figs.  29  and  30  for  192°K  data,  a  correction  of  11  to  16°K  would  be 
necessary  for  the  case  of  ten  rotational  lines  used  in  the  standard  con¬ 
dition  data  and  Fig.  27.  This  yields  rotational  temperatures  from  140 
to  145°K.  For  the  case  where  15  rotational  lines  are  used.  Fig.  23 
would  predict  a  correction  of  19  to  24°K,  and  from  Fig.  27  one  obtains 
values  of  rotational  temperatures  from  138  to  143°K.  For  the  off- 
standard  condition  data,  when  ten  rotational  lines  are  used.  Fig.  23 
predicts  a  correction  of  about  12°K,  which  yields  from  Fig.  28  a  tem¬ 
perature  of  98°K.  When  15  rotational  lines  are  used  Fig.  23  predicts  a 
correction  of  20  to-25°K,  and  this  gives  rotational  temperatures  from 
92  to  97°K.  Again,  reasonable  agreement  between  measured  and  pre¬ 
dicted  Tjj  values  was  obtained  by  taking  into  account  the  effects  of 
density  and  number  of  rotational  lines  used. 

The  (0,  1) /( 1,  2)  band  intensity  ratio  in  the  off-standard  condition 
was  measured  to  be  3.  94  (average  of  ten  integrations)  and  after  correc¬ 
tion  for  wavelength  sensitivity  was  3.  78.  From  Fig.  15  this  gives  a  Tv 
of  1630°K.  An  average  of  another  ten  integrations  gave  a  corrected 
value  of  3.  72  for  the  (0,  1)  /( 1,  2)  band  intensity  ratio,  and  this  gives 
a  Tv  of  1665°K.  An  average  of  eight  integrations  of  the  (0,  2)/(l,  3)  band 
ratio  gave  an  uncorrected  ratio  of  2.  16  and  a  corrected  value  of  2.  09, 
which  gives  from  Fig.  15  a  Tv  of  1740°K.  Therefore,  the  average  value 
of  Tv  for  the  off-standard  condition  is  1678°K  with  a  standard  deviation 
of  ±56°K,  which  is  to  be  compared  with  the  T0  value  of  1662°K  at  which 
temperature  the  specific  heat  ratio,  y,  of  nitrogen  is  1.  31.  Assuming 
one -dimensional  isentropic  flow,  if  the  vibrational  mode  remained  in 
equilibrium  to  the  nozzle  throat  and  remained  frozen  afterwards,  one 
would  find  the  free -stream  vibrational  temperature  to  be  1444°K  which 
when  compared  with  the  measured  value  of  Tv  of  1678  ±  56°K  indicates 
that  the  vibrational  mode  remains  frozen  in  energy  very  close  to  its 
reservoir  value  throughout  the  flow. 

The  average  intensity  of  the  (0,  0)  band  during  flow  at  the  off- 
standard  condition  relative  to  the  intensity  at  0.  100  torr  at  room  tem¬ 
perature  was  measured  to  be  1.54,  which,  when  used  with  the  average 
Tv  =  1678°K  in  Fig.  18,  gives  pm  =  2.  64  x  10"7  ±  0,  06  x  10"  7  gm/cm^. 
This  density  value  is  within  15  percent  of  the  predicted  value  of 
2.  32  x  10-7  gm/cm^;  however,  for  this  condition  it  is  difficult  to  blame 
beam  spreading  for  the  discrepancy. 

In  the  standard  flow  condition  an  average  of  eight  integrations  gave 
a  (0,  1) /( 1 , .2)  band  intensity  ratio  of  2.  91  and  a  corrected  value  of  2.  79, 
which  yields  a  Tv  of  2125°K  from  Fig.  15.  An  average  of  eight 
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integrations  gave  a  (0,  2)/(l,  3)  ratio  of  1.  50  and  a  corrected  value 
of  1.  45,  which  yields  a  Tv  of  2600°K.  Therefore,  the  average  value 
of  Tv  for  the  standard  condition  is, 236 3  ±  237°K,  which  is  to  be  com¬ 
pared  with  the  reservoir  temperature,  T0  =  2365°K,  for  which  y  of 
nitrogen  is  1.296.  For  one -dimensional,  is entropic  equilibrium  flow, 
Tv  at  the  throat  of  the  nozzle  would  be  2060°K  which  would  be  the  free- 
stream  value  of  Tv  if  the  vibrational  mode  froze  at  the  nozzle  throat. 

As  in  the  previous  case,  the  Tv  data  indicate  vibrational  freezing 
throughout  the  flow  which  is  not  unexpected  in  view  of  the  low  reservoir 
pressure.  Time  was  not  available  for  density  determination  in  the 
standard  condition. 

A  spectral  scan  of  the  beam  induced  radiation  in  the  off-standard 
condition  from  3500  to  6000  A  is  shown  in  Fig.  31. 


5.3  10  N2  NOZZLE-TUNNEL  L 

The  10  N2  nozzle  is  a  contoured  nozzle  designed  for  arc -heated 
flow,  the  parameters  of  which  are  given  below: 

M,  -  10 

Standard  Condition 
1.  22  atm 
3100°K  (5580°R) 

145°K  ( 26 1°R) 

2444  m/sec  (8020  ft/sec) 

6.  16  x  10”®  gm/cm3 
(3.  85  x  10'6  lbm/ft3) 
p  =  0.  0205  torr 
Mffl  =  10.  15 
ERTP  =  0.  041  torr 

The  nozzle  may  also  be  operated  in  a  cold-flow  mode  to  achieve  a  con¬ 
dition  of  an  extremely  lowr  rotational  temperature. 

The  diameter  of  the  flow  core  was  approximately  2.  0  in. ,  and  beam 
measurements  were  made  on  the  centerline  1  in.  from  the  nozzle  exit. 
Normally  a  0.  25 -in.  beam  length  was  observed  by  the  spectrometer. 
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It  was  expected  on  the  basis  of  laboratory  temperature  measure¬ 
ments  at  192°K  (Ref.  7),  the  experiments  of  Robbin  and  Talbot  (Ref.  13) 
(Fig.  32),  and  the  predicted  low  free-stream  density  for  this  nozzle 
that  no  more  than  a  10°K  correction- would  be  necessary  to  eliminate 
both  the  density  and  kinetic  temperature  dependencies  of  the  electron 
beam  determined  rotational  temperature.  Therefore,  of  all  the  flow 
fields  available  for  study  in  Tunnels  L  and  M,  it  was  felt  that  the  10  N2 
condition  of  Tunnel  L<  was  the  most  ideal. 

The  results  of  three  sets  of  six  rotational  line  scans  are  shown  in 
Fig.  33,  and  a  typical  spectral  scan  is  shown  in  Fig.  34.  The  variation 
of  Tr  with  the  number  of  rotational  lines  used  was  not  inordinately  large 
when  compared  with  previous  data  and,  in  fact,  is  similar  to  the  192°K 
laboratory  data  shown  in  Fig.  29.  However,  the  allowable  10°K  cor¬ 
rection  to  the  data  yields  a  temperature  from  175  to  180°K,  when  from 
6  to  15  rotational  lines  are  used,  which  is  30  to  35°K  above  the  pre¬ 
dicted  value.  -Approximately  36  more  rotational  scans  yielded  results 
similar  to  those  shown  in  Fig.  33. 

The  average  value  of  the  (0,  1) /( 1 ,  2)  band  intensity  ratio  was  2.  19 
corrected  to  2.  10,  which  gives  a  Tv  of  3280°K.  The  average  (0,  2)/(l,  3) 
band  intensity  ratio  was  1.  215  corrected  to  1.  175,  which  gives  a  Tv  of 
349 0°K.  The  average  value  of  the  (0,  3)/(l,  4)  band  intensity  ratio  was 
0.  8095  corrected  to  0.  830,  which  gives  a  Tv  of  3440°K.  The  average 
value  of  the  (1,  3)/(2,  4)  band  intensity  ratio  was  1.  764  corrected  to 
1.  706,  which  gives  a  Tv  of  3620°K.  Therefore,  the  average  value  of  Tv 
is  3458  ±  141°K  which  is  358°K  higher  than  the  reservoir  temperature, 

T0,  of  3100°K. 

The  average  intensity  of  the  (0,  0)  band  during  flow  relative  to  the 
intensity  at  0.  100  torr  at  room  temperature  was  measured  to  be  0.  314, 
which,  when  used  with  the  average  Tv  =  3458°K  in  Fig.  17,  gives 
pa  =  6.  49  x  10"§  ±  0,6  x  10~8  gm/cm3.  This  is  in  good  agreement  with 
the  predicted  value  of  6.  15  x  10"8  gm/cm^. 

A  scan  of  the  beam  induced  radiation  in  the  standard  condition  from 

O 

3500  to  6000  A  is  shown  in  Fig.  35  which  shows  that  all  significant  radia¬ 
tion  is  attributable  to  excited  or  ionized  nitrogen. 

The  10  N2  nozzle  was  also  operated  in  a  cold-flow  condition.  The 
results  of  two  sets  of  six  rotational  line  scans  are  shown  in  Fig.  36, 
with  a  typical  spectral  scan  shown  in  Fig.  37.  The  variation  of  Tr 
with  K'  is  considerable  above  K'  =  6.  Predicted  values  of  the  flow 
parameters  at  this  condition  are  not  known  at  this  time,  and  definition 
of  the  isentropic  core  by  impact  pressure  surveys  has  not  been  obtained. 
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5.4  18  N2  NOZZLE-TUNNEL  M 

The  18  N2  nozzle  is  contoured  and  designed  for  arc-heated  flow, 
the  parameters  of  which  are  given  below: 

Ma  *  18 

18.  9  to  19.  2  atm 
2800  to  3200°K  (5040  to  5760°R) 

44  to  52°K  (79.2  to  93.6°R) 

2438  to  2621  m/sec  (8000  to  8600  ft/sec) 

5.  45  x  10“8  to  6.  25  x  10“8  gm/cm^ 

( 3.  4  x  10“ 8  to  3.9  x  10" 8  lbm/ft8) 
p,,,  =  0.  060  to  0.  068  torr 
=  17.  9  to  18.  3  , 

ERTP  =  0.  036  to  0.  041  torr 

Beam  measurements  were  made  at  two  axial  positions:  15.  24  cm  (6  in. ) 
and  31.  12  cm  (12.  25  in. )  downstream  of  the  nozzle  exit.  The  diameter 
of  the  flow  core  was  approximately  25.  4  to  30.  5  cm  ( 10  to  12  in.  ),  and 
a  2.  0-cm  beam  length  was  observed. 

Normally,  four  rotational  spectra  scans  were  used  in  each  rota¬ 
tional  temperature  determination.  Figures  38  and  39  show  typical  rota¬ 
tional  temperature  results  on  the  centerline  at  the  12.25-in.  and  6.0-in. 
position,  respectively,  and  Fig.  40  shows  a  typical  rotational  spectral 
scan. 

Unlike  Tunnel  L,  Tunnel  M  has  the  capability  of  on-line  computer 
calculation  of  tunnel  flow  parameters  using  the  measurement  of  impact 
pressure  p^  j  therefore,  for  each  rotational  spectral  scan,  a  tunnel 
parameter  datum  point  was  taken.  Tunnel  M  was  also  operated  at  two 
slightly  different  mass  flow  rates,  but  the  calculated  flow  parameters 
were  always  within  the  range  given  in  the  preceding  list.  The  dif¬ 
ferent  flow  rates  will  be  referred  to  as  miow  and  n^high* 

Figure  41  shows  the  results  of  a  rotational  temperature  radial  pro¬ 
file  measurement  at  the  x'  =  12.  25-in.  position  with  miow.  The  tem¬ 
perature  profile  is  seen  to  be  flat  within  ±5°K  out  to  12  cm.  The  rota¬ 
tional  temperatures  for  this  profile,  as  for  all  Tr  profiles,  were 
determined  using  ten  rotational  lines.  The  average  rotational  tem¬ 
perature  over  the  radial  distance  of  0  to  12  cm  was  68.  4°K,  which  can 
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be  corrected  to  59.4°K  from  Fig.  23  using  a  calculated  average  of 
6.  4  x  10'8  gm/cml  However,  since  the  predicted  free-stream  tem¬ 
perature  was  51.  2°K,  the  rotational  temperature  determined  is  approxi¬ 
mately  16  percent  higher  than  the  calculated  free-stream  temperature. 

Figure  42  shows  the  results  of  a  Tr  radial  profile  measurement  at 
the  x7  =  6.  0-in,  position  with  m^ow  The  temperature  profile  is  ob¬ 
served  to  be  flat  within  ±5°K  in  the  radial  direction  out  to  14  cm.  The 
average  Tr  over  the  14  cm  radial  distance  was  62.  3°K,  which  can  be 
corrected  to  53.  3°K  using  Fig.  23,  which  agrees  to  within  9  percent 
with  the  average  aerodynamic  free-stream  temperature  value  of  49.  1°K. 

Figure  43  shows  the  results  of  a  Tr  radial  profile  at  the  x7  =  6.0-in. 
position  with  mhigh*  The  temperature  profile  is  again  observed  to  be 
flat  out  to  14  cm  within  ±7°K  and  the  average  Tr  over  this  distance  was 
60.  6°K,  which  can  be  corrected  to  51.6°K.  This  is  in  good  agreement 
with  the  average  aerodynamic  Tm  of  49.  6°K. 

An  average  value  of  Tv  was  obtained  at  x'  =  12.  25  in.  on  the  center- 
line  using  the  four  band  intensity  ratios  (0,  1)/(1,  2),  (0,  2)/(l,  3), 

(0,  3) /( 1,  4),  and  (1,  3)/{ 2,  4)  yielding  Tv  =  3216  ±  58°K  (Table  V).  The 
average  TQ  for  which  the  Tv  data  were  obtained  was  2825  ±  50°K;  there¬ 
fore,  a  difference  of  nearly  400°K  exists  between  the  measured  value 
of  Tv  and  the  calculated  value  of  TQ. 

The  Tunnel  M  test  cabin  could  not  be  evacuated  below  0.  75  torr 
with  no  flow,  so  a  density  calibration  at  room  temperature  could  not  be 
made.  However,  the  radial  variation  of  free-stream  density  relative 
to  the  centerline  density  was  measured.  Figure  44  shows  the  radial 
variation  of  density  measured  with  the  electron  beam  and  also  pitot 
probe  at  x7  =  12.  25  in. ,  rh^ow.  It  is  shown  that  the  electron  beam  data 
predicts  a  smaller  flat  density  core  than  the  pitot  probe  data,  but  the 
observed  difference  may  be  attributable  to  the  electron  beam  data 
having  been  obtained  at  a  different  test  cabin  pressure  p,p,  as  is  dis¬ 
cussed  in  Section  6.  1.  2.  Figure  45  shows  the  radial  density  profile  at 
the  x7  =  6.  0-in.  position,  m^0w»  anc*  radia.l  density  gradient  is  ob¬ 
served  to  be  much  less  as  one  moves  closer  to  the  nozzle  exit  for  the 
rh^ow  condition.  Figure  46  shows  the  radial  density  profile  for  the 

x7  =  6.  0-in.,  riikigk  condition.  The  profile  is  observed  to  be  very 
similar  to  the  x7  =  12.  25-in. ,  m^ow  condition  with,  however,  a  some¬ 
what  smaller  flat  core. 
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SECTION  VI 

CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FUTURE  WORK 


6.1  SUMMARY  AND  CONCLUSIONS 

A  summary  of  the  Tunnels  L  and  M  electron  beam  data  is  shown  in 
comparison  with  the  calculated  free -stream  parameters  in  Table  VI. 

6.1.1  Vibrational  Temperature 

From  Table  VI  it  can  be  seen  that  all  values  of  Tv  except  one  meas¬ 
ured  by  the  electron  beam  technique  are  higher  than  the  calorimetrically 
inferred  reservoir  temperatures,  T0,  and  all  values  are  significantly  higher 
than  the  predicted  temperature  at  the  nozzle  throat,  which  Tv  of  N2  should 
approach  if  significant  vibrational  energy  transfer  occurs  in  the  expansion 
from  the  reservoir  to  the  throat,  and  that  the  difference  in  the  experi¬ 
mental  values  of  Tv  and  T0  increase  as  TQ  increases.  There  are  three 
obvious  possible  sources  of  this  discrepancy  which  will  be  discussed. 

6.1. 1.1  Inaccurate  Calorimetric  Determination  of  T0 

If  one  assumes  equilibrium  reservoir  conditions,  i.  e. , 

Tv  (experimental)  =  Tr  =  T0,  it  can  be  seen  that  the  error  in  T0  and 
reservoir  enthalpy  for  all  flow  conditions  except  the  off-standard  9  N2 
nozzle  flow  amounts  to  at  least  10  percent,  which  is  more  than  5  per¬ 
cent  larger  than  the  experimental  uncertainty  associated  with  the 
calorimetric  determination  of  TQ,  as  described  in  Ref.  12.  For  this 
reason  this  source  of  discrepancy  is  considered  no  further. 

6.1. 1.2  Inaccurate  Excitation  Band  Strength 

It  has  been  stated  that  calculable  excitation  Franck-Condon  factors 
q(v'  ,  VJ*  )  have  been  used  in  Eq.  (11)  rather  than  the  presently  unknown 
excitation  band  strengths  p(v',  vj*  ),  and  their  usage  is  justified  on  the 
basis  that  during  the  short -duration  (=  10"  ^  sec)  high  energy  collision 
of  a  primary  electron  with  a  nitrogen  molecule,  the  molecule  has  little 
time  to  move  either  by  rotation  or  vibration  and  the  Franck-Condon 
principle  is  applicable.  However,  if  sufficient  low  energy  secondary 
electrons  exist  in  the  flow,  a  significant  fraction  of  excitations  could 
occur  for  which  the  ratio  of  interaction  time  and  period  of  vibration  is 
considerably  increased;  this  situation  would  violate  the  criterion  of 
applicability  of  the  Franck-Condon  principle  and  would  necessitate  usage 
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of  the  excitation  band  strengths  p(v'  ,  ),  for  which  the  first-order 

approximation  can  be  written  as 

p(v',  Vj")  =  q  (v',  Vi")[l  +b  fv'iVl"]  (14) 

where  b  is  an  unknown  excitation  band  strength  parameter.  If  Eq.  (11) 
is  now  calculated  with  b  as  a  parameter,  the  b  =  0  curve  represents 
the  Franck-Condon  approximation,  and  if  the  source  of  the  Tv  discrep¬ 
ancy  can  be  attributed  to  the  usage  of  the  b  =  0  curve  given  by  Eq.  (11), 
one  can  attempt  to  effectively  scale  the  experimental  values  of  Tv  to  an 
estimated  free -stream  Tv  by  the  judicious  selection  of  a  unique  value 
of  b. 


For  the  low  pQ  expansions,  i.  e. ,  pQ  <  2  atm,  it  is  reasonable  to 
expect  Tv  to  freeze  near  the  reservoir  value  TQ,  and  one  can  deter¬ 
mine  the  value  of  b  required  to  reduce  Tv  (experimental)  to  Tv  *  TQ. 
Doing  this  with  the  (0,  2) /( 1,  3)  band  ratios  for  the  9  N2  off-standard 
condition  one  finds  b  =  0.  75  and  for  the  9  N2  standard  condition  one 
finds  b  =  1.3,  whereas  the  10  N2  data  give  b  =  1.  5.  For  the  18  N2  noz¬ 
zle  of  Tunnel  M,  b  is  found  to  be  approximately  2.  0.  However,  if  one 
assumes  that  Tv  =  T  at  the  throat  of  the  18  N2  nozzle  because  of  the 
high  reservoir  pressures  in  Tunnel  M,  then  one  gets  b  >  3.  Therefore, 
the  values  of  b  thusly  obtained  vary  from  0.  75  to  3,  yielding  no  unique 
value  of  b  to  enable  a  scaling  of  the  electron  beam  results  of  Tv.  This 
result  does  not  imply  that  b  =  0,  or  that  Franck-Condon  excitation  is 
appropriate,  but  rather  that  a  breakdown  of  Franck-Condon  excitation 
rules  is  not  the  sole  source  of  the  discrepancy.  For  all  flow  fields 
investigated  b  is  always  greater  than  or  equal  to  0.  75,  so  that  it  is 
entirely  possible  that  an  increase  of  b  from  0  to  0.  75  is  indicated  and 
that  an  additional  source  (or  sources)  of  discrepancy  exists.  In  the 
event  b  =  0.  75,  the  experimental  Tv  from  the  (0,  2) /( 1,  3)  bands  is 
changed  from  1740  to  1662°K  for  the  off-standard  9  N2  nozzle  where 
T0  =  1662°K.  For  the  standard  9  N2  condition,  if  b  =  0.  75,  Tv  is  de¬ 
creased  from  2600  to  2500°K,  and  the  corresponding  TQ  =  2365°K;  the 
10  N2  result  for  Tv  in  the  test  section  is  changed  from  3490  to  3325°K, 
whereas  T0  =  3100°K  and  the  18  N2  result  for  Tv  is  decreased  from 
3280  to  3100°K,  while  T0  =  2825°K.  Therefore,  it  can  be  seen  that  for 
three  of  the  four  cases,  the  decrease  in  Tv  effected  by  adjusting  b  from 
0  to  0.  75  is  insufficient  to  account  for  the  finding  of  Tv  greater  than  TQ. 
However,  it  should  be  noted  that  each  of  these  three  discrepancies  then 
becomes  less  than  10  percent. 
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6.1. 1.3  Electron  Heating  of  the  Vibrational  Mode  of  N2 

As  a  final  cause  of  the' (presumed)  Tv  anomaly  one  may  question 
the  assumption  of  the  attainment  of  thermal  equilibrium  in  the  stilling 
chamber,  or  reservoir,  thereby  enabling  Tv  to  have  a  higher  starting 
value  for  expansion  than  T0.  The  mechanism  whereby  the  nonequilib¬ 
rium  steady-state  condition  Tv  4  T0  in  the  reservoir  is  maintained 
could  well  be  the  low  energy  electrons  produced  in  the  arc  heater, 
carried  into  the  reservoir  gas  and,  by  means  of,  say,  negative  ion 
formation  and  subsequent  ionization,  producing  an  excessively  high 
vibrational  mode  energy;  the  collision  can  be  represented  as 

(v=0)  +  e“  -»  (N2)-  -*  N2(v  =  1,  2,  3 . )  +  e~ 

For  the  case  of  the  Tunnel  M  18  N 2  flow  one  would  expect  thermaliza- 
tion  to  occur  in  the  reservoir  since  the  residence  time,  tr,  of  an  N2 
molecule  in  the  reservoir  is  over  1000  times  the  magnitude  of  the  vibra¬ 
tional  relaxation  time,  tv,  of  N2  at  3000°K;  one  would  reasonably  expect 
this  excessive  vibrational  energy  to  be  transferred  by  collisions  into 
the  rotational  and  translational  modes  before  emergence  into  the  nozzle. 
The  only  way  this  cannot  be  the  case  is  for  electrons  to  persist  within 
the  reservoir  without  being  thermalized,  thereby  maintaining  a  non¬ 
equilibrium  situation  such  that  Tv  4  T0  in  the  settling  chamber. 

The  next  question  to  answer,  assuming  this  supposition  to  be  true, 
is  whether  calorimetric  determinations  of  T0  are  sufficiently  sensitive 
to  detect  the  nonequilibrium  environment.  If  Tv  in  the  reservoir  were 
4000°K,  maintained  by  electron  heating,  and  both  the  translational  and 
rotational  modes  had  temperatures  of  3000°K,  the  enthalpy  difference 
associated  with  this  assumed  value  of  Tv  and  the  case  where 
Tv  =  TQ  =  3000°K  is  approximately  65,  5  cal/gm,  as  compared  with  the 
932  cal/gm  total  enthalpy  of  the  reservoir;  this  difference  of  approxi¬ 
mately  7  percent,  which  is  marginal  with  regard  to  detection  by  the 
calorimeter,  is  sufficiently  large  to  explain  the  anomalies  in  the  Tv 
measurements. 

Therefore,  of  the  three  sources  of  discrepancy  discussed,  it  is 
concluded  that,  although  no  unique  adjustment  of  b  resolves  the  Tv 
anomaly,  an  increase  of  b  from  0  to  0.  75  coupled  with  the  assumption 
on  a  nonequilibrium  reservoir  condition,  Tv  4  T0,  can  as  easily  account 
for  the  discrepancy  as  can  the  assumption  that  b  =  0  and  that  the  in¬ 
equality  of  Tv  and  T0  is  the  sole  source  of  the  result.  Neither  approach 
invalidates  accuracy  statements  of  Ref.  12  regarding  the  calorimetric 
evaluation  of  reservoir  enthalpy. 
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6.1.2  Density  Measurements 

Free -stream  density  measurements  in  the  Tunnel  L  flows  were 
always  higher  than  the  .predicted  values.  However ,  for  all  cases  the 
measured  values  were  within  15  percent  of  the  predicted  values,  and  it 
is  believed  that  a  large  portion  of  this  error  can  be  attributed  to  in¬ 
accurate  measurement  of  beam  current  caused  by  beam  spreading 
effects.  It  is  to  be  noted  that  such  a  systematic  discrepancy  in  p  might 
be  attributed  to  the  Tv  dependence  of  the  density  measurement,  as 
shown  by  Figs.  17  and  18,  particularly  since  the  experimental  values 
of  Tv  were  consistently  higher  than  expected.  Figure  47  shows  the 
ERTP  as  a  function  of  Tv  for  a  constant  band  intensity  value,  approxi¬ 
mately  equal  to  the  measured  value,  for  the  Tunnel  L  10  N2  nozzle.  To 
lower  the  experimental  ERTP  by  approximately  5  percent  to  the  pre¬ 
dicted  ERTP  requires  a  decrease  in  the  measured  Tv  («  3500°K)  of 
approximately  300°K  which  is  not  unreasonable.  However,  for  the 
Tunnel  L  9  N2  nozzle,  off-standard  condition.  Fig.  48  shows  that  the 
required  decrease  in  p  of  approximately  15  percent  requires  that  the 
free-stream  Tv  be  less  than  1200°K  rather  than  the  measured  value  of 
1678°K,  and  a  change  in  Tv  of  this  magnitude  is  felt  to  be  unwarranted. 
Consequently,  spreading  effects  of  the  electron  beam  remain  as  the 
primary  suspected  cause  of  the  discrepancy  in  density  values. 

The  radial  density  profiles  measured  for  Tunnel  M  consistently 
showed  a  smaller  uniform  density  core  diameter  than  that  shown  by 
impact  pressure  measurements.  Furthermore,  the  much  steeper  pro¬ 
file  at  x'=  6.  0  in.  for  m^ig^  than  for  m^ow  is  quite  interesting,  for  this 

slight  difference  in  m  values  should  have  had  no  effect  on  the  profile. 
However  these  differences  cannot  be  viewed  seriously  because  it  is 
possible  that  the  ratio  of  tank  pressure  to  free-stream  pressure  may 
not  have  been  constant  in  all  cases;  during  the  density  data  acquisition 
for  the  x'=  6.  0  in. ,  m^g^  condition,  a  malfunction  occurred  in  the  data 
system  that  caused  the  loss  of  tunnel  flow  parameter  data. 

6.1.3  Rotational  Temperature 

Rotational  temperature  measurements  for  all  conditions  except  the 
10  N2  nozzle  and  18  N2  nozzle  at  x'=  12.  25  in.  were  within  10  percent 
of  the  predicted  free-stream  temperature  values.  For  the  18  N2  noz¬ 
zle  at  x'=  12.  25  in.  the  measured  value  of  Tr  was  16  percent  higher 
than  the  predicted  value.  However,  if  the  proper  density -temperature 
correction  factor  for  50°K  rotational  temperature  measurement  were 
used  instead  of  the  78°K  correction  factor,-  as  given  in  Fig.  23,  this 
16  percent  error  will  quite  likely  be  reduced  as  will  the  small  discrep¬ 
ancy  detected  at  the  6-in.  position.  Unfortunately,  for  the  35  to  40°K 
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discrepancy  between  measured  and  predicted  rotational  temperature  for 
the  10  N2  nozzle,  there  would  appear  to  be  little  further  correction 
possible  unless  rotational  relaxation  effects  are  surprisingly  large.  The 
discrepancy  is  glaring  and  should  be  of  concern  to  electron  beam  experi¬ 
menters  and  gas  dynamicists. 

In  an  effort  to  examine  the  effects  that  the  electron  beam  measured 
values  would  have  on  the  calculated  flow  parameters  for  the  10  N2  noz¬ 
zle,  D.  E.  Boylan,  ARO,  Inc.,  AEDC,  has  calculated  four  different 
cases  for  10  N2  flow  conditions.  The  standard  flow  parameters  were 
calculated  using  either  standard  (aerodynamic)  tunnel  measurements, 
electron  beam  measurements,  or  combinations  of  the  two  as  known 
values  in  gas  dynamic  equations.  Case  I  (Table  VII)  conditions  were 
calculated  using  the  standard  tunnel  measurements  of  mass  flow  rate, 
impact  pressure,  reservoir  pressure,  total  enthalpy,  and  product  of 
p.U.  obtained  by  a  probe.  Case  II  conditions  were  calculated  assuming 
the  flow  parameters  to  be  defined  by  the  most  plausible  measured  value 
of  rotational  temperature,  the  impact  pressure,  the  reservoir  pressure, 
and  the  nozzle  mass  flow  rate.  To  produce  such  a  flow  would  require  a 
500°K  higher  value  of  T0  and,  moreover,  an  effective  throat  diameter 
of  0.  154  in.  —  significantly  larger  than  the  measured  value  of  0.  148  in. 
Case  HI  conditions  were  calculated  assuming  the  flow  parameters  to  be 
defined  by  the  electron  beam  measured  values  of  rotational  temperature 
and  density  and  the  standard  measurement  of  mass  flow  rate  and  reser¬ 
voir  pressure.  For  Case  III  a  value  of  impact  pressure  approximately 
20  percent  higher  than  that  normally  measured  would  be  required,  and 
it  is  the  belief  of  the  gas  dynamicists  that  this  much  error  in  impact 
pressure  measurement  is  impossible.  For  Case  IV  conditions  the  elec¬ 
tron  beam  measured  value  of  density  and  the  standard  measurements  of 
impact  pressure,  reservoir  pressure,  and  mass  flow  rate  were  assumed 
to  define  the  flow  parameters.  Although  this  approach  does  not  eliminate 
the  gap  between  measured  rotational  temperature  and  predicted  free- 
stream  temperature,  the  other  quantities  tabulated  do  not  differ  by  im¬ 
possible  amounts  when  the  Case  I  and  Case  IV  results  are  compared. 

Note  also  that  in  Case  II  a  result  is  obtained  for  TQ  which  differs 
markedly  from  the  value  derived  by  directly  measuring  HQ,  and  in 
Case  III  a  serious  discrepancy  arises  when  the  computed  p£  is  compared 
to  the  directly  measured  value.  Thus,  it  may  be  said  that  Cases  II 
and  III  are  not  nearly  as  well  supported  by  impact  pressure  and  reser¬ 
voir  enthalpy  data  as  is  Case  IV,  but  for  Case  IV  the  calculated  value 
of  T^  differs  significantly  from  the  measured  value  of  rotational  tem¬ 
perature. 

Using  a  different  rotational  line  intensity  data  reduction  technique 
called  the  least-squares  quadratic  fitting  technique,  all  the  rotational 
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spectral  data  were  reanalyzed.  Results  from  this  method  of  data  reduc¬ 
tion  are  also  shown  in  Table  VI.  Good  agreement  is  obtained  with  the 

* 

experimental  values  obtained  from  the  least -squares  linear  fitting  tech¬ 
nique.  The  quadratic  fitting  technique  is  discussed  in  detail  in  Appen¬ 
dix  III. 


6.2  RECOMMENDATIONS  FOR  FUTURE  WORK 

For  the  10  N2  nozzle  the  discrepancy  in  measured  and  predicted 
temperature  should  be  resolved.  Two  experiments  should  be  performed: 
(1)  various  stilling  chamber  lengths  should  be  used  to  see  if  the  Tr 
discrepancy  is  a  result  of  nonequilibrium  processes  in  the  stilling  cham¬ 
ber,  and  (2)  the  10  N2  free-stream  conditions  should  be  duplicated  with 
a  conical  nozzle  operating  at  a  much  lower  stagnation  temperature.  This 
latter  experiment  was  brought  to  mind  since  the  departure  of  Tr  from 
the  predicted  value  has  been  observed  to  increase  with  Tq  as  shown  in 
Figs.  49  and  50.  This  variation  has  also  been  observed  by  Sebacher 
(Ref.  14)  and  MacDermott  and  Marshall  (Ref.  15).  However,  MacDermott 
and  Marshall  were  able  to  obtain  agreement  between  rotational  tempera¬ 
ture  measurements  and  static  temperature  predicted  using  finite-rate 
theory  with  magnified  vibrational  rates  based  on  their  vibrational  tem¬ 
perature  measurements.  Such  an  approach  is,  of  course,  impossible 
for  the  10  N2  nozzle  data  since  measured  vibrational  temperatures  were 
always  greater  than  or  equal  to  the  predicted  T0. 

Future  flow  diagnostic  efforts  in  Tunnels  L  and  M  should  also  be 
concentrated  on  direct  measurements  of  T0  and  Tv  of  the  reservoir 
using: 

1.  Spectral  line  reversal  apparatus 

2.  Emission  spectroscopic  observations  of  the  N2+  [l“] 
system  for  Tv 

3.  Consideration  of  seeding  N2  with  trace  amounts  of  H2 
or  He  to  look  for  Stark  broadening  of  emission  lines 
attributable  to  high  values  of  the  electron  density  in 
the  reservoir, 

4.  Seeding  the  flow  with  trace  amounts  of  CO  or  CN  to 
determine  the  reservoir  values  of  Tr  (or  T0)  using 
molecular  band  emission  attributable  to  arc  heater 
excitation. 

5.  Investigation  of  the  feasibility  of  making  spatially  re¬ 
solved  measurements  of  density  and  temperature  in  the 
reservoir  by  means  of  laser  Raman  scattering. 
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A  repeat  of  all  the  pa  measurements  using  the  electron  beam  at  a 
higher  energy  would  .also  be  desirable  in  order  to  investigate  the  effect 
on  the  measurements  reported  herein. 
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Fig.  1  Elevation  View  of  Tunnel  L 
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Fig.  2  Elevation  View  of  Tunnel  M 
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Fig.  4  Schematic  Drawing  of  Electron  Beam  Vacuum  System 
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Fig.  5  Schematic  Drawing  of  Tunnel  L  Electron  Beam  Installation 


Fig.  6  Schematic  Drawing  of  Tunnel  M  Electron  Beam  Installation 
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Fig.  7  Diagram  of  Beam  Receiver  Cup  and  Electrical  Connections 
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Fig.  8  Tunnel  M  Auxiliary  Beam  Receiver  Cup 
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View  Port  (Glass) 


Fig.  9  Schematic  Drawing  of  Tunnel  L  1-Meter  Spectrometer  Setup 
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Electron  Gun  Housing 


Fig.  10  Schematic  Drawing  of 


M  1  -Meter  Spectrometer  Setup 
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Fig.  12  Tunnel  M  Electronic  Integrator 
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Relative  Sensitivity 


O  Tunnel  L  Spectral  Sensitivity 
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Fig.  13  Spectral  Sensitivity  of  Spectrometer  Optical  System  for  Tunnel  L 
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Relative  Sensitivity 
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Fig.  14  Spectral  Sensitivity  of  Spectrometer  Optical  System  for  Tunnel  M 
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Fig.  15  Theoretical  Vibrational  Band  Intensity  Ratios  versus  Vibrational  Temperature 
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Fig.  16  Theoretical  Vibrational  Band  Intensity  Ratio  versus  Vibrational  Temperature 
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Relative  Intensity,  (o,  i)  Band 
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q(v vf ')  Calculated  by  Hornkohl 
Vibrational  Level  Population  from  Muntz  (Ref.  3) 

Tv.  °K 


Nj  Pressure,  torr  at  298°K 

Fig.  17  (0,  i)  Band  Relative  Intensity  as  a  Function  of  Density  (ERTP)  and  Tv 
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Fig.  19  Typical  (0,  0)  Band  R-Branch  Rotational  Line  Spectral 
Scan,  4  ISI2  -Nozzle,  Standard  Condition 
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o  4/28/69,  Rotational  Spectral  Scan  DWL  0112  -  0117,  4  N2  Nozzle 
A  4/28/69,  Rotational  Spectral  Scan  DWL  0118  -  0123,  4  N2  Nozzle 


0  2  4  6  8  10  12  14  16  18  20  22 
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Fig.  20  Rotational  Temperature  versus  Number  of  Rotational  Lines  Used  in  the  Temperature  Determination, 
4  N2  Nozzle,  Standard  Condition 


o  5/8/69,  Rotational  Spectral  Scan  DWL  0124  -  0129,  4  N2  Nozzle 
a  5/8/69,  Rotational  Spectral  Scan  DWL  0130  -  0135,  4  N2  Nozzle 
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Fig.  21  Rotational  Temperature  versus  Number  of  Rotational  Lines  Used  in  the  Temperature  Determination, 

4  N2  Nozzle,  Off-Standard  Condition 
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R' 'actual 


°  4  N2  Nozzle,  0. 270  p0  Condition 

Q  Extrapolated  Laboratory  Data,  Tactua|  =  78°K,  p  ~  3. 6x  10"^  gm/cm3  (Ref.  7) 

0  Laboratory  Data,  Tactua|  ’  78°K,  P  a  1. 8  x  10"  ^  gm/cm3  (Ref.  7) 

v  Laboratory  Data,  Tactua| s  78°K,  P  =  3. 0  x  10"8  gm/cm3  (Ref.  7) 


Fig.  22  Ratio  of  Measured  Rotational  Temperature  to  Actual  or  Predicted  Temperature  versus  Number 
of  Rotational  Lines  Used  in  the  Temperature  Determination,  4  N2  Nozzle  Data  and  78°  K 
Laboratory  Data 
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Fig.  23  Empirical  Prediction  of  Rotational  Temperature  Measurement  as  a  Function  of  Nitrogen  Density 
and  Number  of  Rotational  Lines  Used  for  Temperatures  Near  78°  K 
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Fig.  24  Electron  Beam  Induced  Radiation,  4  N2  Nozzle,  Standard  Condition 
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Fig.  25  Typical  (0,  0)  Band  R-Branch  Rotational  Line  Spectral  Scan,  9  l\l2  Nozzle, 
Standard  Condition 


Fig.  26  Typical  (0,  0)  Band  R-Branch  Rotational  Line  Spectral  Scan,  9  N2  Nozzle 
Off-Standard  Condition 
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o  5/29/69,  Rotational  Spectral  Scan  DWL  0148  -  0153,  9  N2  Nozzle,  Standard  Condition 
*  5/29/69,  Rotational  Spectral  Scan  DWL  0154  -  0159,  9  N2  Nozzle,  Standard  Condition 


Fig.  27  Rotational  Temperature  versus  Number  of  Rotational  Lines  Used  in  the  Temperature 
Determination,  9  N2  Nozzle,  Standard  Condition 
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o  5/12/69,  Rotational  Spectral  Scan  DWL  0136  -  0141,  9  N2  Nozzle,  Off-Standard  Condition 
a  5/12/69,  Rotational  Spectral  Scan  DWL  0142  -  0147,  9  N2  Nozzle,  Off-Standard  Condition 
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Fig.  28  Rotational  Temperature  versus  Number  of  Rotational  Lines  Used  in  the  Temperature 
Determination,  9  N2  Nozzle,  Off-Standard  Condition 


o  N2+(l"),  (0,  0),  Pc  -  0.0044  torr, 


K'(l  through  K'  Spectral  Lines  Used) 


Fig.  29  Ratio  of  Measured  Rotational  Temperature  to  Chamber  Wall  Temperature  versus  Number  of  Spectral  Lines 
Used  for  the  Temperature  Determination  with  Chamber  Pressure  as  a  Parameter 
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Fig.  30  Empirical  Prediction  of  Rotational  Temperature  Measurement  as  a  Function 
of  Nitrogen  Density  for  Temperatures  Near  192°  K 
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Fig.  31  Electron  Beam  Induced  Radiation,  9  N2  Nozzle,  Off-Standard  Condition 
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Fig.  32  Ratio  of  Actual  Temperature  to  Measured  Rotational  Temperature 

versus  Actual  Temperature  for  Densities  on  the  Order  of  10®  gm/cm3 
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a  Rotational  Spectral  Scan  DWL  0028  -  0033,  10  N2  Nozzle,  Standard  Condition 

o  Rotational  Spectral  Scan  DWL  0045  -  0050,  10  N2  Nozzle,  Standard  Condition 

a  Rotational  Spectral  Scan  DWL  0064  -  0069,  10  N2  Nozzle,  Standard  Condition 
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Fig.  33  Rotational  Temperature  versus  Number  of  Rotational  Lines  Used  in  the  Temperature  Determination, 

10  N2  Nozzle,  Standard  Condition 
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Fig.  35  Electron  Beam  Induced 
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ion,  10  N2  Nozzle,  Standard  Condition 
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Fig.  37  Typical  (0,  0)  Band  R-Branch  Rotational  Line  Spectral  Scan,  10  N2  Nozzle, 
Cold  Blow  Condition 
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o  3/2/70,  Rotational  Spectral  Scan  DWM  0001  -  0006,  C.  L.,  (0,0)  Band,  x'=  12.25  in. 

•  4/13/70,  Rotational  Spectral  Scan  DWM  0037  -  0041,  C.  L.,  (0,  0)  Band,  x'-  12.25  in. 
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Fig.  38  Typical  Rotational  Temperature  Plots  versus  Number  of  Rotational  Lines 
Used  in  the  Temperature  Determination,  18  N2  Nozzle,  x'  =  12.25  Inches 
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o  5/8/70,  Rotational  Spectral  Scan  DWM  0089,  C.  L.,  (0,  0)  Band,  x'=6.0  in. 

•  6/4/70,  Rotational  Spectral  Scan  DWM  0116  -  0119,  C.  L,  (0,  0)  Band,  x'=  6.0  in. 
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Fig.  39  Typical  Rotational  Temperature  Plots  versus  Number  of  Rotational  Lines  Used  in  the  Temperature 
Determination,  18  N2  Nozzle,  x’  =  6.0  Inches 
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Fig.  40  Typical  (0,  0)  Band  R-Branch  Rotational  Line  Spectral  Scan,  18  N2  NOzzle 
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Fig.  44  Relative  Nitrogen  Density  Radial  Profile,  18  N2  Nozzle,  x'  =  12.25  Inches,  rh|OW 
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Fig.  45  Relative  Nitrogen  Density  Radial  Profile,  18  N2  Nozzle,  x'  =  6.0  Inches,  m,ow 
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Fig.  46  Relative  Nitrogen  Density  Radial  Profile,  18  N2  Nozzle,  x'  =  6.0  Inches,  mhigh 
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TABLE  I 

PROGRAM  TROT 


- PRt)6Ri\M  r  i-j  m 

TYPF  NU.  NUNIjRrt,  NSUMX Y  ,  iMSUMASO 

TYPE  I MTtbEH  VU.VL 

Ui*EN*>IUN  I  <22>  *G(22>  «  X  Y  t  <j2 )  »  AS'i  I  22 )  tUtLlAY  (  22  >  »  V  f<  22  >  » 

**C^Z);PfHl2Zl,«AWEii2).  lift"  2S)  »  SIGmaTa < lV*  2U>>  . 

_ *_  _  _  tm < i oo ) »  sig(1uo)»  s ( i o n >  •  ssaoojj_jr uooi» 

1  YY(lco) » '"‘z"(ioo>  •  <mnui*  wY'Uotj) »"  lo(ift)” ' 

3o 1  read  33^t  t*i».  lily _ 

IV  (EOF * b  )  331*332 

_311_STi^P _ 1  _ 

332  LUI.U® 

DU  334  L=’<2S _  _ 

Tw ( i  )  =  Twuf 

334  CONTINUE _ 

M  =  0 

1  RtAQ  2f)Q.kUN.  VU.VL _ 

LULU®  LULIJ  ♦  1 

N=  ij 

M  =  M  +  1 
PHInT  loo 
PHInT  l."l 
PHInT  .1  12 

PRINT  lu3. HUN, VU.VL 

PRINT  1  n4 _ _ 

Ht AO  201 *KMAX 

KKK  =  KMfl X _ 

DO  3  K=1»KMAX 

3  HEAP  2o2, I  (K) _ 

TWv,=2oo« 

4  SUMX®.  . 

SU-Y®-. 

SUwXY*-; . _ * 

SU*iXSO«o. 

BU=2.i)H3“0.02QMVU*.5) _ _ 

BL®1.932“O«020*<VL*.5) 

DO  5  K=l  ,  KMAX _ _ _ 

NU=2.9979?5E*10*  <2546 l .5  ♦  241 9 . 84*'<  Vlji ,5> -23.  19*  <  < VU*.5>**2> 

»-22n7.19» (VL+.5)  +lk  .  1  4«  (  (  VL*  .  S>  »  »2  >  ♦  K«  (K*  1 )  «bU-K«  I K- 1 )  »t>D _ 

NUNOHm=nU/< 2. 9979?5F*1 0*25566. 0> 

WAVE  (K)=2«997925E*1  o/NU*  1  .  pb.*  0  8 _  _ _ 

G<K)  =  <K*EXPF  (5.  7?5B*K/THG>  ♦ <K+i )  »EXPF  <-5*72:>«*  (K*  D  >f  HbJ  >  /  (  2«*k*  i ) 

Y(K)=LOGF(I (K)/< <NUN0HM»«4)«G(K)«*) )  _ ; _ ' _ 

X<K)=1.*K*(K*] > 

XY<k)=X(K)«Y<K) _ 

XSQ(K>=X(K)*XU> 

SUMX=SUMX»X (K) _ 

SUMY=SUMY+Y (K) 

SUMXYgSUMXY*XY(K) _ . _ _ 


5  SUMXSu*SUMXSQ*XSQ (K) 
SOSUHX=SUmX«SUwX 


PHODSUM=SUMX*SUMY 

nsumxy*kmax*sumxy 

NSUmX-SQ=KMAX«SUMXSU 

SLOPES InSUMXY-PRodSUM) / (NSUMXSQ-SUSUMX) 

OROCEP®  <  SiJMXSQ*SUMY-SUMX*SUMXY )  /  <  NSUMX5Q-SQSUMX  > 

TR=- 2.86 29/SLOPE 
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TABLE  I  (Continued) 


SUHj)ELSrj= 

DO  ft  K® It KMAX 

Y  T  (K)  =SLUPE*X  cKTTuHiicr  P 

UELTA.Y(K)=Y(K)-YT(K) _ 

PtW  (K)  =  1  I  AY  (K ) / YT  (K)' 

ft  SUMi)ELSiJ=SUM|)ELSu*Ij(;LTAYCK>*OELTaY  (K_)  _ _ 

SiftHAY=lSuMUtLSU/CKwiAX-2.i  )*».5 

SIbMAr=(2.86?9*SlGMAY/ (SLOPE*SLOPd) ) * ( (KMAX/ (NSUMXSQ-SuSUmX) >**.5) 
TtbTTH=AriSF (TH-Trg) /TR 

IF  (TESTlR.LT.  H.7 _ 

I  Tki,=  TR 

bO  TO _  _ _ 

B  N=N*1 

IF  (N.GT.|)  11.9 _ 

V  Ou  10  K  =  It  KMAX 

1 1,  PRINT  lob.Kt  I  (K) ,G(K) tY (K) ,X(K)  .WAVE  (K)  .PE8(K> 

PRINT  l.-b 

PRINT  ] Q  7 .PR(1I)SUm,NSUmXY  ,  SuSUmX.NSUmXSQ  _ 

PRINT  1  h.slope.orucep 

print  1  ,9 _ _ 

PRINT  110, KMAX, fH.SlGMAT 

_ [RR  (MtK  -QAx)  =  TR _ _  _ 

S1'iMATA(M,  KMAX)  =  SI8MAT 

_ KMAX«KMaX-] _ 

(4U  TO". 

n  print  hi.kmax. ir.sigm«t _ 

TuR ( M  t  K  M  A  X )  =  I R 

SlbMATA (M.  KMAX)  =  SIGMAT _  _ _ 

KMAX=KMAX-'l 

I F  (KmaX.gT.S )  A. 400 _ 

401)’  IF  (LULU. LT. LILY)  1,12 

12  CONTINUE _ 

KM  X  —  KKK 

RRInT  ft  i  _ _ _ _ 

6S  FORMAT’  ( l h  1 1 

IS)  CONTINUE _ _ _ 

IF  (KmaX  .GT.  ST  IB,  300 

18  COnTIndF. _ 

PRINT  T 

_ UU  J  3  J_  =  1  ,  _M_  _ _ _ 

PRINT  H  t  THR< jtKMAX)  ,  TW<J)'t  SIGMA  FA  ( Jt  KMA*>‘ 

_i_3  CONII.mUF  _ _  _  _ 

oo  i4  j  =  it  m 

Y  ( J)  =  IHH  ( Jt  KM  jO/TwU) _ 

S(J)  =  SIf;MrtTA(Jt  KMAX)/TW(  J) 

Sb(J>  =  S(J)**2» 

Y  V  (  j )  =  Y  ( J)  /Sb  (.))  '  . . 

Z(J)  =  t./SSU) 

jJ  =  'Ss(J)/S'S(.j) 

14  »  Y  ( j)  =  t)  ( J)  o  W(J)  •  SS(J) _ 

SUM l  =  0.0  s  SUM?  :  0.0  )  SUM3  =  0.0  !>SUm4  =  0.0 
DO  lb  J  =  It  M 

SUM'V  ="  su«‘l"t  VY(  J) 

SUN?  a  bUM?  ♦  /(J) _ _ 

SUM  1  =  ♦  if Y  ( j) 

lb  bU-i»  =  SUt<‘4  ♦  W  <  j ) _ 
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TABLE  I  (Concluded) 


TRBaR  =  SUMi/SUMg 

S1GBAR  =  (SUM3/(SUM4  *  SUM4 ) )  **  0.5 
PRINT  30*  THBAH*  SlGRAR 

PRINT  210,  KMAX*  M _ 

KMAX  =  KMAX  “  I 

PRINT  pj _ _  _  _ 

23  FORMAT  (IhO*  777) 

GO  TO  1 9 _ 

210  FORMAT  (1H0*  5X,  5HUSING,  13*  18H  LINES*  AVtRAGE  OF*  12*  bh  RUNS> 

30  FORMAT  ( 1H0 * 15HAVERAGE  TEMP  =  *  F15.9*  5X*  jbHAVERAGE  SIGMA  -  * _ 

1F1  i.5) 

70  FORMAT  (1H  *  2()HmEASURED  TEMPERATURE*  SX*  16HWALL  TEMPERATURE, 
*ltiX,  1AHSTANOARO  DEVIATION  //) 

Bfl  FORMAT  (1H  *  3X*  F15.9*  7X*  F15.9*  13X,  Fin. 5) _ 

100  FORMAT  (1h1»/»47X*25HP0TaTi0NAL  BOLTZMANN  PLllT) 

101  FORMAT  (1H  */*57X*5HUSlN0)  _ 

1  o2  FORMAT  <1H  ,/,52X»  14H-UJNTZ  6-FACTOR) 

1 03  FORMAT  (lH  *//* 2l X* 1 2HKUN  NUMBER  =* IX* AB*5X* 1 1 , 1H* *  1 1 ,5H-dAND) _ 

1{>4  FORMAT  (lH.  »//,6X,1hK,  7X*  4Hi  (K)  ,  )  lX»  4HG  (K)  ,  lbX»4HY  (K  >  *  15X  ,  4riX  (K)  * 

«13X*1vHWAvELENQTh*9A*?7H(OELTA  Y(K))/YT(K),  PERCENT) _ 

l.i5  FORMAT  (1H  »<*X*  l2»6X,15,8X,Ell.4»BX,Eli  .4»3X,E11.4,8X,E11.4» 

»l2X*Ell  ,4) _ _ 

1;:6  FORMAT  (1H  ,//,bX,nHSUM(X)*SUM(Y),6X,9HN*SUM(XY),8X, 

»13HSUMCX) aSUM(X) *6X*1  .HN»SUM tX»X) )  _ _ _ 

1  .i7  FORMAT  (1H  ,5X,Ell.4*GX,Ell.4,7X,Eil.4,7X,E11.4) 

1,'lB  FORMAT  (1H  *//,2lX*7HSL0RE  =  .  1 X*  E 1 1  .4  *  1  OX  *  1  lHlNTERCEPT  g*lX*El1,4) 
1  f|9  FORMAT  (lh  .//) 

110  FORMAT  (1H  *.15X*  I2*gH-LINES*5X,24HWOTaT1QNAL  TEMPERATURE  =* _ 

*E 1 1 .4*5X*26HSTANDAH0  DEVIATION  OF  TR  =,Ell.4> 

111  FORMAT  (1H  *15X*I2*6H-LINES*5X,24H _ zl _ 

♦E11.4*5X*26H  "  =,E11.4) 

200  FORMAT  (Aa*lX*Il*lX*Ii) _ 

201  FORMAT  ( I 2 ) 

202  FORMAT  (I«5> _ 

333  FORMAT  (Flfl.3*  12) 

enu  trot _ 
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TABLE  II 

VIBRATIONAL  BAND  STRENGTHS  FOR  EMISSION 


0 

1 

2 

3 

4 

0 

0.6348 

0. 2706 

0.  0784 

0.0192 

0. 0042 

1 

0. 2580 

0.  2179 

0.  2983 

0.  1472 

0. 0505 

2 

0.  0319 

0.  3452 

0.  0501 

0.  2372 

0.  1836 

TABLE  III 

FRANC K-CONDON  FACTORS  FOR  EMISSION 


BB 

0 

1 

2 

— 

3 

4 

0 

0.  6570 

0. 2546 

0. 0686 

0.0157 

0. 0329 

1 

0.  2982 

0.2316 

0.  2838 

0. 1299 

0. 0418 

2 

0. 0482 

0.4065 

0.  0561 

0.  2284 

0.  1635 

TABLE  IV 

FRANCK-CONDON  FACTORS  FOR  EXCITATION 


0 

1 

2 

0 

0. 8794 

0. 1064 

0.  0127 

1 

0. 1183 

0.6876 

0.  1606 

2 

0.  0022 

0. 2012 

0.  5587 

84 


AEDC-TR-71-61 


TABLE  V 

SUMMARY  OF  TUNNEL  M  VIBRATIONAL  BAND  INTENSITY  DATA 


1 

Band  Pair 

Raw  Intensity- 
Ratio 

Corrected  Intensity 
Ratio 

Tv,  °K 

(0,  1)1(1,  2) 

2.  27 

2.  16 

3140 

(0,  2) /( 1,  3) 

1.  25  . 

1.  22 

3280 

| 

(0,  3)/(l„4) 

0. 84 

0.  87 

3221 

(1,  3)/( 2,  4) 

•  » 

1.  77 

1.91 

3222 
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TABLE  VI 

SUMMARY  OF  TUNNEL  L  AND  TUNNEL  M  ELECTRON  BEAM  DATA 


Nozzle 
Flow  Field 


4  N2  Nozzle 
Standard  Conditions 


4  N2  Nozzle 

Off-Slaiidard 

Conditions 


9  N2  Nozzle 
Standard  Conditions 


0  N2  Nozzle 

Oil-Standard 

Conditions 


10  N2  Nozzle 
Standard  Conditions 


10  N2  Nozzle 
Cold  Flow 


18  N2  Nozzle 
x  =  12.  25  in. 


18  N2  Nozzle 
x  6.  0  in. 


"calc. 

°K 


67 


66 


135 


93 


145 


51.2 


49.  1  to  49.6 


^Rmcas.’ 

°K 


60  to  67 
10  Lines 
64  to  68 
15  Lines 

60  to  67 
10  Lines 
64  to  68 
15  I  .ines 


140  to  145 
10  Lines 
138  to  143 
15  Lines 


08 

10  Lines 
92  to  97 
15  lines 


175  to  18 0 
6  to  15  Lines 


59.4 


51.  6  to  53.  3 


^Rmeas.  Q.  fit* 
"K 


66.4 


66. 1 


1 4  J .  7 


03.4 


180.! 


30 


59.  1 


51.  2  to  54.  3 


Tv 


meas. 

K 


2363 


1678 


3158 


3216 


To- 

°K 


298 


298 


2365 


1662 


3100 


298 

2800  to  3200 
2825 


2800  to  3200 


™meas.’ 

gm/cm^ 


p">'  calc.  ’ 
gm/cm^ 


5. 61  x  10 


-7 


4.  8  x  1()’7 


6.56  x  10 


-7 


5.  6  x  10 


-7 


2.02  x  10 


-7 


2.  64  x  10 


-7 


6.49  x  10* 


2.  32  x  10 


6. 16  x  10 


-7 


5.45  x  10'°  to 
6.25  x  10'8 


5.  45  x  10~8  to 


6.  25  x  10 
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TABLE  VII 

TUNNEL  L  10  N2  FLOW  CONDITIONS  FOR  FOUR  CASEST 


Parameter 

Units 

Case  I 

Case  II 

T* 

°K 

2656 

3121 

2845 

2448 

p* 

psia 

9.  756 

9.  770 

9. 761 

9.  749 

T  ) 

0  real 

i 

°K 

3050  . 

(H0  known) 

3573 

3265 

2813 

10.  10 

10.  10 

9.  71 

10.  10 

To, 

°K 

149 

175 

known 

175 

known 

137 

P„ 

M  Hg 

21.  1 

21.  1 

26.  8 

21.  0 

lbm/ft^ 

3. 968  x  10-6 

3. 379  x  10-6 

4.  1  x  10-6 
known 

4.  1  x  lO-6 
known 

u. 

ft/sec 

8248 

8939 

8594 

7909 

Re^ 

per  inch 

392 

317 

389 

437 

K. 

in. 

0. 0389 

0.  0479 

0. 0377 

0.  0347 

% 

psf 

4.  195 

4.  196 

4.  924 

4.  170 

p " 

*o 

MHg 

2780 

known 

2780 

known 

3337 

2780 

known 

d* 

in. 

0.  148 
known 

0.  154 

0.  151 

0.  145 

P  u 

lbm/ft  2-sec 

0. 03273 
known 

0.  03020 

0. 0369 

0.03393 

m 

lbm/hr 

7.  76 

known 

7.  76 
known 

7.  76 
known 

7.  76 
known 

po 

psia 

18 

known 

18 

known 

18 

known 

18 

known 

Case  Ir  Standard  Measurements  Used 

Case  II:  Electron  Beam  Results  (T,,,)  and  p^,  m,  and  pQ 

Case  III:  Electron  Beam  Results  and  p0  and  m 

Case  IV:  Electron  Beam  Results  (p^)  and  pQ',  m,  and  pQ 


^  All  four  calculations  assume  the  flow  is  in  equilibrium  in  the  stilling 
chamber  and  becomes  frozen  at  the  sonic  point  which  is  assumed  to  be  d*. 
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TABLE  VIII 

PROGRAM  TROTQUAD 


PROGRAM  TROTQLAD 

TYPE  INTEGER  VU.  Vl.  CODE 

TYPE  nniJRLF  XL*  Y*  C.  VARY.  VC 

DIMENSION  PHUOIi  HP?).  K(22>.  CODEIA).  *N(22>,  6(2?)  ,  XL  <22. A) 

*  C(6).  VC(6.ft).  X (?2) .  Y(?2) 

1  READ  2000.  (PH(N) ,N=1,10) 

read  2ml .  run.  lines,  vo.  vl 
DO  3  Net. LINES 

re  An  ?m2.  I  in)  .  k  (M 

IF  (K(N).FO.n)  2,3 
?  K (N) =  N 

3  rONTINUE 

nLL  WAVE  (WN.  K.  VU.  VL*  LINES) 

4  PEAD  2003.  (CODE (N) »N=1 .6) 

IF  (EOF. 5'))  S.6 

5  5TOP 

6  TF  (CDOEd  >  .FQ.01  1,7 

7  THGa  200«P 

?  nO  P  N®1  ,1.  INFS 

G<N)  =  (  K(N)  *  FXPF (5,7258  *  K ( N)  /  TpG  )  ♦ 

*  (  K(N)*1  >  •  ExpE (  -S.7?5«  *  (K(N)+i)  /  TRG  >  >  /  (  2*K(N)*1  ) 

Y  (  N )  a  l.((GE(  I  (I'D  /  (  f!(N)*K  (N)  *I»N  (N )  **«  )  ) 

O  X(N)=  KIN)  *  (K (M) * » ) 

CALL  SFTUP  (XL,  MATRIX,  X,  LINES,  CnnE) 

CALL  PEGRFSS  (C,  Vary.  VC,  XL,  Y,  MAtdIX,  LINES) 

TR=  -?.S629  /  C ( ?) 

TFSTTR=  A«SF(  (Tp-TwG)/T»  ) 

IF  (Tf.STTM.GT.1. PE-04)  10*11 

10  TRiis  TR 
GO  TO  « 

11  PRINT  1000,  (PH(M) ,MS),10) 

PRINT  1001,  RUN,  VU,  VL 
PRINT  1002 

JMAX®  MATRIX  -  1 

no  14  Nal.LINFS 

YCAL*  C(l) 

no  ]3  JsI.JMAX 

IF  (CODE ( J* 1 )  .NE  .<))  1?,13 

12  YCAL=  YCAL  ♦  C ( J* 1 )*X(N)**j 

13  CONTINUE 

PER*  100.0  *  (  (Y(N)-YCA|.)/Y(N)  ) 

14  PRINT  1003,  K(N),  I(N),  X ( N) ,  YIN),  YrAL,  PER 
PRINT  inOA,  (CODE (N) ,,=1 .fi) 

PRINT  1005,  VARY 
PRINT  100* 

PRINT  1007,  (C(N> ,N=1 .MATRIX) 

PRINT  100B 

PRINT  1007,  (VClM.N) ,N=1, MATRIX) 

SOTRS  P.BA29  *  Vc(2,2)  /  (  C(2)*C(2>  ) 

PRINT  100P,  TR,  SOTR 
GO  TO  A 

.1660  FORMAT  (IHl,  I0X,  1 0 An) 

1001  FORMAT  (IHfl,  20X,  6HR')N  a  ,  AR,  SX,  If,  IN,,  II,  SH-BAND) 

1002  FORMAT  (1H0,  RX.  1 HK ,  4X.  | HI .  7X,  1HX,  9X,  1 HY ,  ]]X,  4HYCAL* 

*  7X*  AHYEHROR) 

1003  FORMAT  ( 1H  ,  7X,  12,  ?X,  I5,  2X,  F6.1,  2X,  Ell. 4,  2X,  Ell. 4, 

*  ) 

1Q0*  format  mho*  10*«  612) 

1005  FORMAT  (IHfl,  10X,  2EHCTANDARD  OFVTATlON  IN  Y  a,  Ell. 4) 

1006  FORMAT  ( 1  ho  ,  loX,  20HCOEFF I  Cl  ENTS  IN  YCAL) 

1007  FORMAT  (1H  ,  1  OX ,  5(1X.E1I.4>> 

1008  FORMAT  (IHfl,  10X,  3rHSTANOaRO  DFVIaTIonS  In  COEFFICIENTS) 

1009  FORMAT  ( 1 H 0 ,  10X,  4HTR  =.  F11.4.  5X,  ^((STANDARD  DEVIATION  IN  TP 

.  *  *  Ell, 4) 

2000  FORMAT  (  1 0 AR  ) 

2001  FORMAT  (AR,  IX,  Tp,  IX,  II,  IX,  ID 

2002  FORMAT  (IS,  IX,.  TP) 

2003  FORMAT  (  M2  ) 
end  Trot 
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TABLE  VIII  (Continued) 


S'triRDlIT  T  N-  HAVF  (HIM,  K  •  Vll»  V  lL  • 
ntMfNSinN  v»M(?2).  k<?2)*  rViJfioi 
TTRF  IMTFRKK  VII,  V| 

rvii(u=  2./,o'>9^o?t  o-»  s  cvii(?)=- 
rVi>(3)=-l.?2-iS9335E  oi  *  rV0(4>  = 
rvi)<5)=-5.i2t>ul27ftE-ni  i  rvu<ft>= 
rV'l(7)=-2.?0<i93070E-0o  *  CVU(fl)  = 

rVil(9)*-1.0lu)C44?E-OA  4  nVU<llJ)s 

/VLI1I*  2.207247loE  O'*  *  rVL(?l=- 

rvi  (31  =-?,4741  T4R9E-n?  4  CVL<4)=- 
ORIMIJ*  2.0H«511OftE  f»n  $  rRU(?)=- 
f«'J(3)=  3.70‘i2633ftE-0O  *  CHU(4)=- 

r«i'(S)=  6.177S1017E-0*  t  rRU(ft)*- 
rPU{7)=  l.x93ai730t-0» 

0HL<1>=  l.R31926ft3E  on  *  CHL(?>=- 
CVI.I3IS  3.272S155RE-04  4  r«L(4)=- 

rp|.(5)=  4,47731 1  *n£-o? 

Oil*  2S441.50 

XU*  1 .  n*vn  ♦  <1.3 
XT=  l.c 
no  i  Ms) (]n 
XT*  Xr  *  XII 

1  filjs  Gil  ♦  CVIMN)  *  XT 


line  ft) 

.  CV|  (4)  »  CRUI7)  ,  cRL<5> 

-4.4nft54445E  00 
3»3'>R3S594E  00 
4.41S09795E-02 
0.44444393E-03 
•  6.774464S6E-09 
■1 •67009353E  01 
■9,OnM0391E_04 
■3.fHH4ft?A9E-02 
■B.fto44hl39E-04 
■l»7O33«fiS0E-0h 

■2.0oift719bE-02 

■P.aaonBBSlE-OS 


OL=  0.0 

XL  *  1.0*VL  ♦  0.5 
xr*  l.o 
no  2  H:li4 
X  T=  XT  *  X| 

?  GL*  Gl  ♦  CVL(N)  *  XT 
HVU=  CPiK  i  J 
XT*  l.o 
nO  3  M=P,7 
XT*  XT  *  XII 

3  HVU=  HVU  ♦  C«U(Nl  *  XT 
OVL*  CHI.  ( 1  J 

XT*  1.0 

no  4  N=?,G 

XT*  XT  *  X|_ 

4  RVL*  HVL  ♦  CHUM  *  XT 
no  5  n=i.LINFS 

5  WN(M)*  GIJ  *  MVU*K (N)»(K(NJ*1)  -  GL  -  BVL»« ( Nl * ( K ( N ) - 1 ) 

*NORM=  *N(1) 

nO  ft  N=1 .LINFS 
ft  wN(N)*  vIN ( N )  /  XMORm 
RF  TURN 
FMI)  WAVF 
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TABLE  VIII  (Continued) 


SUBROUTINE  SETUP  (XL.  MATRIX.  X,  LINfs.  CODE) 

TYPE  INTEGER  COrvE.  TFMCODF 
TYPE  DOUBLE  XL 

0 1  MENS  I  ON1  XL  (22. M  ,  x(22).  COOElf.).  TFMC00EI6) 
no  1  1*1*22 
00  1  J*1»E 

1  XL ( 1 . J) =  0.0 
MATRIX*  0 

no  i  n=i  «*> 

IF  (CODE  IN) .EQ.l)  2,3 

2  MATRIX*  MATRIX  *  1 

3  CONTINUE 

no  8  1*1 .L1NF5 
no  4  N*1 *6 

*  TEMCOIIE  <N)  *  COOE(N) 

DO  8  J*1*MATRIX 
t)0  7  Ms  l  ,  A 

IF  (TFMCODE(N) .Eo.ll  m.7 

5  XT*  1.0 
kmax*  n-i 

no  «  k*i,kmax 

6  XT*  XT  *  X  ( I ) 

TEMCOOE (N> »  o 
00  TO  8 

7  CONTINUE 

8  XL  ( I*  J)  =  XT 
RFTIJRN 

FNn  SETUP 

SUBROUTINE  REGRESS  (C.  VARY,  VC*  XL.  Y.  MATRIX,  LINES) 
TYPE  DOUBLE  XL.  Y.  A.  Z.  C*  VARY.  VC 

DIMENSION  XL (22, 6> ,  Y<22>.  A(6,6),  Z(*>.  C<6).  VC(6*6> 
DO  2  1*1. MATRIX 
no  2  J*I. matrix 
AlI.JI*  0.0 
no  1  Ks) .LINES 

1  A 1 1 • J) *  A ( I . J)  ♦  XL(K.I)  *  XL(K.J) 

2  A(J.I)*  A(I.J) 

DO  3  1*1. MATRIX 
Z(I)«  0.0 

00  3  J=l. LINES 

3  Z<I)*  Z(D  *  X|.U,I)  *  T ( J) 

CALL  HATINVRT  (A,  MATRIX) 

no  4  1*1. MATRIX 
C(l>*  0.0 
no  4  j*i .matrix 

4  C(I)»  C(l)  ♦  AlI.JI  *  ZtJ> 

vary*  0.0 

DO  5  1*1. LINES 

5  VARY*  VARY  ♦  Y ( I 1  «  YII) 

DO. 6  1*1. MATRIX 

6  VARY*  VARY  -  C<I>  •  Z(I> 

VARY*  A8SF<  VARY/<LINFS-MATRIX)  ) 

no  e  ui.matrix 

DO  7  .J»1»MATRIX 

7  VCd.J)*  A  ( I ,  J )  •  VARY 

B  VC(I.I)*  (  ABSF (  VC(I.I)  )  )**0.5 
RETURN 
ENO  REGRESS 
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TABLE  VIII  (Concluded) 


SUBPOUTINF  MAT1NVRT  (A*  NN) 

TYPF  DOUBLE  A*  AM AX «  PIVOT*  F 
DIMENSION  A  (6*6)  <  LOCATE  (6*3) 
nO  1  N*l*NN 

1  LOCATE (N*  3) a  0 
nO  i 4  N=1«NN 
AMAX=  0.0 

00  6  I=1«NN 

IF  {  LOCATF(1*3).EU.O  I  2*6 

2  00  5  J=1«NN 

•  IF  (  LOCATE! J, 3) .FQ.O  >  3*5 

3  TF  (  ARSF<A(1*J) ) .GT.AMAX  )  4*5 

4  A«AX=  ARSF ! A ( I • J) ) 

TROW*  I 

JC0L=  J 

5  CONTINUE 

6  CONTINUE 

IF  (  AMAX.GT.l'.OF-lS  >  7»  1 8 

7  LOCATE(N.l)*  IROw 
LOCATE  (N*?)«  JCOl. 

LOCATE ( JCOL  *3) *  1 

IF  C  IROW.NE.JC0L  )  8.10 
H  nO  9  J=1 « NN 
SWAPS  A ( I  ROW  « J) 

A( IROW.Jl*  A(JCOL.J) 

9  A (JCOL* J) *  SWAP 

10  PIVOT*  A  (JCOl.*  JCOL) 

A!  JCOL. JCOL)*  l.fl 
no  11  J*1*NN 

11  A  ( JCOL..J)  *  A  ( JCOL*  J)  /  PIVOT 
DO  14  I*1*NN 

TF  (  I. ME. JCOL  I  12*14 

12  F=  A(I.  JCOL) 

A ( I , JCOL) =  0.0 
no  13  J*1*NN 

13  A  ( I  *  J)  *  A(I.J)  -  F*A  (  ICOL*  J) 

14  CONTINUE 

no  17  Nsl.NN 
I  *  nn»n»i 

IF  (  LOCATF (L . 1 ) .NE. LOCATE (L*2)  )  15*17 

15  TPOW*  LOCATE (L* l ) 

•  ICOL=  LOCATE  (L*  2) 
nO  16  K  a 1  *  NN 
SWAP*  A(K.TROW) 

A (K , I ROW) *  A ( K* JCOL) 
a (K * jcol) s  Swap 

16  CONTINIIF 

17  CONTINUE 
RETURN 

18  PRINT  1000 
PFTURM 

i 5 6 o  format  ( i>-  .  iox,  ishsingular  matrix* 

FNO  MAT1NVRT 
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TABLE'  IX 
PROGRAM  TRADV 


PROGRAM  TRADV  . 

nIMFNSION  TR ( 100) .  TW(100>.  SIG  C 1 00 )  »  S(lQO)»  SS<100>»  Y(100>* 

1  YY(100)*  Z<100>.  W(100)»  WY(foO).  iDao) 

PRINT  40 _ 

READ  60»  ID 

PRINT  90 ♦  ID _ _ 

3  REAO  lOtXLINES 

TF(EOF,50)  50,51 _ 

50  STOP 

51  continue _ 

READ  1 0* IA1 

REAP  20,  (TR ( I)  .  TH(1)<  SlG<I>.  I*l*Ifti> 

PRINT  70 

PRINT  BO.  (TR  ( I )  »  Tui(T)  «  SlG(I)  ,  I=l»iAl) 

10  FORMAT  (13) 

20  FORMAT  (2F15.9.  FI 0 .5) 

30  FORMAT  ( 1H0* 15HAV£RAGF  TEMP  =  ,  F15.9,  5X.  16HAVERAGE  SIGMA  =  ♦ 

1F10.5) _ 

40  FORMAT  (1H1 ) 

60  FORMAT  (10A8) _ 

70  FORMAT  <1H  ,  20HMEASU»ED  TfMPERATURE » "sX.  16HWALL  TEMPERATURE, 
«TqX,  IRHSTANDARD  DEVIATION  //) 

80  FORMAT  <1H  ,  3X.  F15.9.  7X,  F15.9*  13x,  FlO. 5) 

90  FORMAT  (25 X,  1QA8 ) _ 

00  l  J=  1  ♦  I A1 

YU)  s  TR(J)/TW(.J)  _  _ 

S(J)  =  SIG(J)/TW (J) 

_ SS(J)  =  S(J) **2. _ 

YY(J)  =  Y(J)7sS(.J) . . 

_ 7U)  =  1./SSIJI 

W(J)  =  SS(i)/SS(J> 

1  W Y ( J )  =  W(J)  »  W(J)  *  SS(J) _  _  _ 

SUM1  =  0.0 

SUM2  =^0.0 _ _ _ 

SUM3  =  0*0 

_ SUM4  =  Q.O _ _ 

nO  2  K=l , T A 1  . 

SUM  j  a  SUMj  »  YY(K) _ 

SUM2  =  SUM?  ♦  Z(K)  . 

SUM3  =  SUM3  ♦  W_Y(K> _ 

2  SUM4  =  SUM 4  ♦~W(K) 

TRBAR  =  SUM1/SUM? 

SI'GBAR  =  '< SUM3/  ( SUM 4  *  SUM4 )")"  **  0-5 

PRINT  30,  TRBAR.  SlGBflR _ 

PRINT  in.  XLINES 

PRINT  ip,  I A 1 _ _ 

GO  TO  3 
FNP 
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APPENDIX  III 

ROTATIONAL  TEMPERATURE  CALCULATIONS  USING  A  LEAST-SQUARES 
QUADRATIC  FITTING  PROCEDURE 


Ashkenas  (Ref.  8)  and  Williams  (Ref.  7)  have  observed  that  as  the 
number  of  spectral  lines  used  to  determine  Tr  is  reduced  the  deduced 
Tr  agrees  increasingly  better  with  the  true  gas  temperature.  For 
instance,  using  a  typical  rotational  spectral  scan  obtained  in  nitrogen 
at  a  pressure  of  0.  005  torr  at  liquid  nitrogen  temperature  (78°K),  a 
linear  fit  through  13  rotational  lines  yields  a  Tr  =  88.  2°K  with  c r  =  1.2°K, 
whereas  a  similar  fit  through  only  the  first  eight  rotational  lines  yields 
Tr  =  82.  7°K  and  tr  =  1.  5°K.  An  inspection  of  the  Boltzmann  plot  of  the 
intensity  values  shows  this  variation  of  Tr  with  the  number  of  spectral 
lines  to  be  the  result  of  curvature  in  the  dependence  of  y  on  x,  where 
y  and  x  are  as  defined  in  Eqs.  (5)  and  (6).  For  this  case  it  now  appears 
that  y  is  at  least  a  quadratic  function  of  x  instead  of  the  linear  relation¬ 
ship  as  predicted  by  the  Muntz  model,  and  this  dependence  can  be 
written  as 

y  =  B0  +  ajx  +  a2x2  ( III  -  1 ) 

If  one  assumes  that  the  causative  agent  for  this  nonlinearity  does 
not  affect  the  linear  term  in  x,  one  sees,  that  aj,  and  thus  Tr,  is  given 
by  the  first  derivative  of  y  with  respect  to  x,  the  slope  of  y,  evaluated 
at  x  equal  to  zero,  or  K  =  0,  where  it  is  assumed  that  aj  =  hcB0/kTR 
as  before.  It  is  emphasized  that  no  a  priori  justification  can  be  given 
for  such  an  assumption.  An  iterative,  least -squares  quadratic  fit  of  the 
spectral  data  to  Eq.  (Ill - 1)  was  performed  using  Program  TROTQUAD 
presented  as  Table  VIII. 

A  least -squares  quadratic  fit  through  the  13  points  of  the  preceding 
example  yields  Tr  =  77.  6°K  and  cr  =  2.  3°K,  which  is  in  excellent  agree¬ 
ment  with  the  true  gas  temperature. 

For  a  given  spectral  scan  a  rotational  temperature  and  standard 
deviation  are  obtained  from  Program  TROTQUAD,  and  for  a  number  of 
spectral  scans  at  the  same  condition  a  weighted  average  rotational  tem¬ 
perature  and  standard  deviation  are  calculated  using  Program  TRADV 
shown  in  Table  IX.  This  program  uses  Eqs.  (7)  through  (10)  which  were 
explained  previously. 

Applying  this  procedure  to  the  laboratory  data  of  Ref.  8,  good  agree¬ 
ment  between  average  rotational  temperature  and  test  chamber  wall 
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temperature  is  found  provided  that  the  nitrogen  number  density -electron 
beam  current  product,  Nxig  (ma/cm^)  is  greater  than  approximately 

1015  ma/cm^.  Since  all  the  free-stream  conditions  of  Tunnels  L  and  M 
which  were  investigated  had  number  densities  such  that  the  product 
N-piB*  where  ig  is  the  measured  electron  beam  current  used  for  the 
measurement,  was  greater  than  10^  ma/cm^  and  less  than  10^ma/cm^, 
it  was  considered  worthwhile  to  analyze  the  tunnel  data  using  the  em¬ 
pirical  quadratic  fitting  method,  and  the  results  are  shown  in  Table  VI 
where  it  is  seen  that  excellent  agreement  is  obtained  between  the  quad¬ 
ratic  fitting  results  and  the  theoretically  predicted  values  except  for 
the  10  N2  nozzle.  For  the  10  N2  nozzle  condition  a  total  of  54  rota¬ 
tional  scans  were  analyzed  and  averaged.  An  average  rotational  tem¬ 
perature  of  180.  5°K  was  determined,  still  approximately  35°K  above  the 
predicted  value. 

Using  both  laboratory  and  tunnel  data  the  coefficient  a2  (Eq.  ( III  - 1 ) ) 
is  plotted  in  Fig.  Ill  - 1  as  a  function  of  temperature,  and  it  is  seen  that 
a 2  is  a  strong  function  of  temperature  and  can  be  inferred  to  be  only 
weakly  dependent  on  N2  density  over  this  temperature  range,  assuming 
NxiB  ^  10^5  ma/crn^. 

It  can,  therefore,  be  concluded  that  for  nitrogen  over  the  tempera¬ 
ture,  density,  and  beam  current  ranges  investigated,  the  quadratic 
form  of  y  as  a  function  of  x  is  a  more  accurate  representation  than  the 
linear  for  all  cases  except  the  10  N2  flow,  and,  additionally,  this  func¬ 
tional  relationship  yields  values  of  aj,  the  functional  form  of  which  is 
hcB0/hTp^,  which  are  independent  of  N2  density.  However,  caution 
must  be  exercised  for  applications  to  any  data  lying  outside  the  investi¬ 
gated  regions  and  for  which  the  gaseous  composition  is  other  than 
pure  N2. 
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